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We use cluster Dynamical Mean Field Theory to study the simplest models of correlated elec- 
trons, the Hubbard model and the t-J model. We use a plaquette embedded in a medium as a 
reference frame to compute and interpret the physical properties of these models. We study various 
observables such as electronic lifetimes, one electron spectra, optical conductivities, superconduct- 
ing stiffness, and the spin response in both the normal and the superconducting state in terms of 
correlation functions of the embedded cluster. We find that the shortest electron lifetime occurs 
near optimal doping where the superconducting critical temperature is maximal. A second critical 
doping connected to the change of topology of the Fermi surface is also identified. The mean field 
theory provides a simple physical picture of three doping regimes, the underdoped, the overdoped 
and the optimally doped regime in terms of the physics of the quantum plaquette impurity model. 
We compare the plaquette Dynamical Mean Field Theory results with earlier resonating valence 
bond mean field theories, noting the improved description of the momentum space anisotropy of 
the normal state properties and the doping dependence of the coefficient of the linear temperature 
dependence of the superfluid density in the superconducting state. 

PACS numbers: 71.27.+a,71.30.+h 



I. INTRODUCTION 

The origin and the nature of superconductivity in 
strongly correlated materials is one of the greatest chal- 
lenges in modern condensed matter theory. It received 
renewed attention with the discovery of the high temper- 
ature superconductivity in copper oxide based materials. 
While these materials have been studied intensively over 
the past decades there is still no concensus as to what 
are the essential physical ingredients responsible for the 
high temperature superconductivity phenomena and how 
it should be modelecH>Md > 5 i 6 l 7 i 8 l 9,io 1 ii 1 i2 l i3 1 ^, 

P.W. Anderson proposed that the high temperature 
superconductivity phenomena was intimately connected 
to the proximity to a parent Mott insulating stat e 15 : 16 . 

Developing precise connections between the proxim- 
ity to a Mott insulator and high temperature supercon- 
ductivity has proved to be a difficult problem. Sug- 
gestive conclusions have been reached using slave boson 
method a 17 : 18 , variational wave function a 19 ' 20 , and gauge 
theory techniques^. However, lack of theoretical tools has 
made difficult to prove that simple models are sufficient 
to explain the phenomenas surrounding cuprates. For ex- 
ample it is still strongly debated whether the existence of 
superconductivity with a high critical temperature and 
a pscudogap is a genuine property of the models studies, 
or, an artifact of the approximations employed to solve 
the model. 

Over the past decade, significant progress in the field 
of correlated electrons has been achieved through the de- 
velopment of Dynamical Mean Field Theor} ^ 21 : 22 . In its 
single site version, this method describes lattice models, 
in terms of a single site impurity problem embedded in 
a medium. The method has been very successful in de- 
scribing, and even predicting numerous properties of a 



large number of materials^^^^^^ 2 *^^. Cluster 
extensions of this method, Cluster Dynamical Mean Field 
Theory (CDMFT) (for reviews seeHlll), have been pro- 
posed and are currently a subject of intensive investiga- 
tions. 

In this article we apply the cluster dynamical mean 
field approach to construct a mean field theory of the 
simplest models of strongly correlated materials, the one 
band Hubbard and t-J models, using a 2 x 2 cluster, 
namely the plaquette as the basic mean field reference 
frame. 

There are several motivations for constructing a mean 
field theory based on a plaquette embedded in a dynam- 
ical bath of conduction electrons : a) A plaquette em- 
bedded in a a self consistent medium can describe the 
physics of singlet formation, which is very important in 
the t-J and Hubbard model. There are two roads of sin- 
glet formation, the Kondo effect, in which a spin can 
form a singlet with a bath of conduction electrons, or, 
the superexchange mechanism which locks two spins on 
a bond in a singlet state, b) A plaquette in a medium 
is a minimal unit to describe d-wave superconductivity 
and antiferromagnetism on the same footing, given that 
their order parameters (as well as that of other forms of 
order competing with superconductivity) naturally fit on 
a plaquette. 

From a methodological perspective, mean field theory 
allows to study physical properties of different phases as a 
function of control parameters, whether they arc stable or 
metastable. For example, we will study the evolution of 
the superconducting state, together with the underlying 
normal state, which appears as a metastable phase be- 
low Tc- From a theoretical perspective metastable states 
are only defined within a mean field theory, but they are 
of clear physical relevance. Furthermore, comparison re- 
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sponsc or correlation functions in both the normal and 
the superconducting state gives important clues as to the 
mechanism of superconductivity. 

A clear understanding of the evolution of well defined 
mean field phases of simplified model is an important 
step towards constructing the phase diagram of realis- 
tic Hamiltonians. Even if a phase is not realized as the 
thermodynamically stable phase in a mean field treat- 
ment of a simplified Hamiltonian, it could be stabilized by 
adding additional longer range terms in the Hamiltonian 
without significantly altering the short distance proper- 
ties described by the mean field theory. Furthermore a 
good understanding of the different mean field states can 
be useful in elucidating the results of numerical studies 
in larger finite clusters, since complicated patterns in a 
finite size system, may be a reflection of phase separation 
among different competing mean field phases. 

The study of minimal models such as the t-J model 
or the Hubbard model, describing a system near a Mott 
transition is an important first step towards understand- 
ing real materials. From a study of minimal models, one 
can learn what aspects follow from just the proximity 
to a doping driven Mott transition. This is a necessary 
step before the importance of other physical effects such 
as the disorder or the electron phonon interactions cer- 
tainly present in the real materials, can be ascertained. 
A basic question yet to be elucidated, is to which extent a 
minimal model of the doping driven Mott transition such 
as the t-J model, describes at the qualitatively level, the 
physical properties of the cuprates. If indeed the qual- 
itative low energy physics of the cuprates results from 
the proximity to a Mott insulating state, as described 
by a minimal model of this phenomena, then the results 
can be refined by including more realistic band structure, 
for example nearest and next nearest neighbor hoppings, 
longer range interactions, disorder and coupling to the 
lattice, as well as incorporating a multiband situation 
which is needed to describe the physics in a wider energy 
range. It is possible to carry out these studies in the 
more realistic framework the combination of electronic 
structure methods with DMFT, a subject which is left 
for future studies. One should also ascertain the size of 
the corrections to the mean field theory, by either ex- 
panding around mean field theory 3 - 3 - or by increasing the 
cluster size 3 - 4 -. 

Several studies have already shown that the Hub- 
bard model treated within cluster DMFT on a 2 
by 2 plaquette successfully describes many proper- 
ties of the high temperature superconductors. For 
example the competition of antiferromagnetism and 
superconductivit y 35 ' 36 ' 37 ' 38 ' 39 , the existence of a pseudo- 
gap at low dopin g 40 ' 41 ' 42 ' 43 ' 44 ' 45 ' 46 , and the formation of 
Fermi arcs 4 ^* 4 ^. 

These phenomena involve short range non local corre- 
lations. In CDMFT the approach to the Mott insulator 
is characterized by the growth of the non local compo- 
nents of the self energy, which is responsible for the phe- 
nomena of momentum space differentiation and the for- 



mation of lines of zeros in the Green's function at zero 
temperature. Surprising manifestations of strong correla- 
tions such as the transfer of optical spectral weight upon 
condensation^ 9 -, the existence of an avoided quantum crit- 
ical point£°- underlying the superconducting dome, and 
the presence of two distinct gap a 51 ' 52 in the supercon- 
ducting state of the undcrdoped cuprates. The approach 
describes well an anomalous incoherent normal stat o 49 ' 53 
which is lifted by the onset of superconductivit y 50 ' 54 . 

Other studies of the Hubbard model using large clus- 
ters at values of U < 8t have focused on the convergence 
of the critical temperature 3 ^. In a series of publications 
it has been shown that the d wave superconducting state 
is well described by spin fluctuation theor y 34 ' 55 ' 56 . To 
which extent the physics of well defined quasiparticles 
interacting with spin fluctuations responsible for the pair- 
ing, can be carried over to strong coupling regime, is an 
important open problem, which can be only be addressed 
by gaining a better understanding of the large U limit of 
the Hubbard model, which is the focus of this article. 

Hence we focus on understanding the physical content 
of the plaquette mean field theory in the regime where 
the interaction strength is large enough to drive a Mott 
transition at half filling with a substantial Mott Hubbard 
gap. We gain insights by comparing the superconducting 
state with the underlying normal state. For example, we 
study the evolution of the Fermi arcs with temperature, 
and trace the mechanism of superconductivity to the op- 
timization of the supcrcxchangc energy. We connect the 
maximum critical temperature with anomalies at opti- 
mal doping, resulting from a maximum in the inelastic 
scattering rate. The techniques introduced in this paper, 
provides a simple interpretation of the cuprate phase di- 
agram in terms of the occupations of a small number of 
cluster eigenstates or pseudoparticles which describe a 
mean-field coarse-grained version of the important exci- 
tations of the lattice system, and we use them to describe 
different experimental probes, tunneling optics, neutron 
scattering, in both the normal and the d-wave supercon- 
ducting phase. The superconducting state is character- 
ized by two energy scales: one increases with decreasing 
doping, and one decreases with decreasing doping. The 
first can be identified with the photoemission gap in the 
antinodal region while the second can be identified with 
the slope of the dirac cone along the Fermi surface. We 
investigate the effect of the latter scale on the penetration 
depth. 

The organization of the paper is the following: In 
section [H] we summarize the formalism and introduce 
the models, the cluster schemes and the impurity 
solvers, i.e. the continuous time Quantum Monte Carlo 
(CTQMC) 57 ' 58 and a generalization of the non crossing 
approximation (NCA^ 9 .^^^. SectionQTQdescribes the 
evolution of the cluster Green's functions and the self en- 
ergies as a function of doping. We identify the existence 
of an anomalous scattering rate describing the nodal re- 
gion of the lattice model, which peaks at a characteris- 
tic doping <52 C in the normal phase. The scattering rate 
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is drastically reduced in the superconducting state. We 
identify a second characteristic doping 8\ c at which an- 
other self energy diverges, and connect this phenomena 
to the formation of lines of zeroes in the Greens function. 

One can view CDMFT in the superconducting phase 
as a generalization of the Migdal Eliashberg theory to 
strongly correlated electron systems, and we present the 
frequency dependence of the superconducting order pa- 
rameter in section IIVI A great strength of the mean 
field theory, is that it allows us to study the "normal" 
state underlying the superconducting state and how it 
evolves with temperature. This is done for the tunneling 
density of states in section IIVI for the optical conductiv- 
ity in [V] and for the magnetic properties in section IVII 
This comparison between the mean field normal state and 
the mean field superconducting state, establishes the su- 
perexchange as the main pairing mechanism, as surmised 
in the RVB theory. 

The pseudoparticles representing plaquette eigen- 
states, are not only technical tools to set up strong cou- 
pling impurity solvers but provide a physical picture of 
the excitations of the system and we use them to interpret 
the CDMFT results in section IVTTl We conclude with the 
connection between our method and an earlier simpler 
mean field theory approach based on the plaquette, the 
slave boson mean field theory and closely related meth- 
ods. For related work advancing the RVB concepts u sing 
single site DMFT on multiorbital models see Refs. I63ll64l 

II. FORMALISM 

In this section we summarize the methodology used 
for our investigation. Two minimal models of the prox- 
imity to a Mott transition were considered: the Hubbard 
model and the t-J model. There are several different 
versions of Dynamical Mean Field Theory For example 
in addition to standard DMFT, an extended version of 
DMFT (EDMFT)^^^ 1 .^^^^^ which replaces all 
the non local terms in the interaction (namely the ki- 
netic energy and the superexchange ) by a fermionic and 
a bosonic bath has been proposed. There are also nu- 
merous variants of cluster Dynamical Mean Field Theory 
which differ by the dynamical medium surrounding the 
plaquette (hybridization function of the impurity model). 
Finally, the solution of the impurity model that results 
from the CDMFT mapping, can be carried out with dif- 
ferent impurity solvers. In this work we use two comple- 
mentary solvers, the non crossing approximation (NCA) 
and a continuous time Quantum Monte Carlo (CTQMC) 
method. 

The goal of this article, is to highlight physical proper- 
ties which follow generally from the proximity to a Mott 
insulating state which are captured by a local approach, 
namely cluster DMFT. For this reason, we have focused 
on physics which emerges from both the Hubbard and 
t-J model, and which is captured by all the different 
cluster schemes (Cellular DMFT—, Dynamical Cluster 



Approximation^! and their extended versions). While 
we mention some quantitative differences between these 
schemes, the stress is on qualitative main conclusions 
that can be obtained with all quantum cluster schemes. 
In order to keep the presentation clear and the article 
relatively concise we provide only methodological details 
which are not available in the literature. To avoid unnec- 
essary duplication, results for a given physical quantity 
are presented with only one cluster scheme and impu- 
rity solver, chosen to demonstrate more clearly a physical 
point. 

A. Models 

One of the more studied models in the field of strongly 
correlated electrons is the Hubbard model defined by the 
Hamiltonian 

ija i 

It consist of a hopping term and an on-site repulsion. 
To be above the Mott transition we take an on-site re- 
pulsion U = 12t. 

A second model of great interest is the t-J Hamiltonian 

" »..<W,' ' ^»<< S < S .<- ( 2 ) 
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It contains two terms; first describes the kinetic energy 
which dclocalizes the holes introduced by doping, and the 
second represents spin-spin interaction. In this work we 
take J/t = 0.3. 

In the t-J model, a constraint forbidding all double oc- 
cupancy must be enforced, and will be treated exactly in 
this work. In the spirit of understanding general features 
of the proximity to the Mott state we include only the 
nearest neighbor hopping t = 1 (t! = 0). 

B. Extended and standard DMFT 

In DMFT, the non local terms in the Hamiltonian cou- 
pling are replaced by a a coupling to a bath of conduc- 
tion electrons. In the Hubbard model the only non local 
term is the kinetic energy, and this leads to the standard 
DMFT mapping which is described in many reviews^. 
In the t-J model also the superexchange interaction con- 
nects different sites, and applying the DMFT philosophy 
to that term also, leads to the Extended DMFT equa- 
tions. 

Here we outline the derivation of the Extended version 
of the cluster DMFT^^ 1 .^^ . We first employ Hubbard- 
Stratonovich transformation to decouple the non-local in- 
teraction term of the t-J model leading to the following 
action 
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Here $ is the Hubbard-Stratonovich vector bosonic field 
which decouples the spin-spin interaction, 

The many body theory described by the action above 
can be summarized in a functional: 

T[g,V] = -TrlogtG^-^-Tr^E] 

+ iTrlog^ 1 - n) + ^Tr[Pn] + $[0,2?] (4) 

Here functional <f> of the exact Baym-Kadanoff func- 
tional contains all two particle irreducible diagrams of 
an electron-boson system with propagators Q, T>. Ex- 
tremizing the functional Eq. (fj| leads to the exact 
Dyson equations for this system. Cluster approxima- 
tions are obtained by restricting the functional to a sub- 
set of trial Greens functions. In the Cellular-DMFT (C- 
DMFT) 23 i 70 , the $ functional is approximated as follows: 
The full lattice is covered by non-overlapping clusters. 
The functional within each cluster is treated exactly, i.e., 
if two lattice points i and j are inside the same cluster 
iC-DMFTjg..^..] = <]»-'••" c;, ; . /•»,/. If however, i and 

j are in different clusters, <J> functional is set to zero. In 
this way, short range correlations within the cluster are 
treated exactly, while long range correlations are ignored. 

Cluster approximations arc obtained by replacing 
the exact functional <f> in Eq. (j4]) by its cluster 
counterpart. The saddle point equations then be- 
come ^cluster = cluster )/SG cluster i ^-cluster — 
— 25^ (D c i us t e r)/ SHduster- 

The fluctuating bosons 3> q in the Extended DMFT 
formalism allow to keep some out-of cluster short range 
correlations and describe better the spin fluctuations by 
allowing the cluster spin to relax more efficiently through 
its direct exchange interaction with the bath. We will see 
that this leads to higher superconducting critical tem- 
peratures. Apart this quantitative difference, we did not 
find any qualitative difference between extended version 
(which employs bosons to describe spin fluctuations be- 
tween the clusters) and the results of the non-extended 
version of DMFT. 



C. Cluster Schemes and impurity models 

There are several cluster schemes in use, in the study 
of correlated electron materials. The Dynamical Cluster 



Approximation (DCA)21 can be thought of as a coarse 
graining in momentum space, obtained by relaxing the 
conservation of momentum. Rather than treating the 
infinite number of lattice k points and corresponding 
Green's functions Qk, the $ functional is approximated 
to depend only on the Green's function of a few clus- 
ter momenta which we will denote by capital letters K 
and Q. The cluster Green's functions of the approximate 
functional ^[^k, 2?q] arc obtained by course graining the 
exact Green's functions, i.e., C?k — * Gk. = SkeK ^ k an< ^ 
£>q — ► T>Q = SqeQ^Q where the sum J2keK is over 
those k momenta in Brillouin zone which correspond to 
certain cluster momenta K (see Ref. HH and Ref. FzTh . 

The results of this paper were obtained with both DCA 
and C-DMFT. Again, all the qualitative features to be 
discussed in the next sections can be seen with both 
methods. Since DCA is a cluster method with a simple 
interpretation in momentum space while C-DMFT has 
a simple interpretation in real space, the fact that the 
qualitative physics emerges from both approaches, sug- 
gests that the physical properties that we discuss in this 
article, are genuine properties of cluster dynamical mean 
field theory on a plaquette irrespectively of the specific 
cluster scheme used. 

We summarize the abbreviations used in the remaining 
of the text: 

• CDMFT: cluster DMFT. 

• C-DMFT: Cellular DMFT-Z2. 

• DCA: Dynamical cluster approximation 7 ^. 

• EC-DMFT: Extended version of Cellular DMFT 

• EDCA: Extended version of Dynamical cluster 
approximatio n 71 ! 72 

A great advantage of all cluster DMFT formulations is 
that the complicated functional Dyson equations for the 
self energies and cluster response functions can be witten 
in terms of an impurity model 
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which is numerically tractable and where the effective 
Weiss fields A and Xo" 1 have to obey the following self- 
consistency conditions 

G = ^(iw-tfk-EO))- 1 = (iuj-E lmp -i:(iuj)-A(iijj)- 1 

k 

(6) 

X = ^(MH + Jq)- 1 = (MM + Xo 1 ^))" 1 (7) 

k 

which merely express the fact that the cluster quanti- 
ties, computed from the impurity model l/(iuj — Ei mp — 
Ti(ioj) — A(iu)), have to coincide with the lattice lo- 
cal quantities when summing over the reduced Brillouin 
zone. Namely, in the C-DMFT, the lattice was divided 
into non-overlapping clusters, hence the summations over 
k run over the reduced Brillouin zone. Here M plays the 
role of the spin self-energy which is computed from the 
local susceptibility and Weiss field by M = x _1 — Xo > 
as evident from the Eq. {7| . 

A special feature of the 2x2 plaquette is worth stress- 
ing: the cluster momentum K is a good quantum number 
and therefore local quantities like Green's function G or 
hybridization A take a diagonal form: 
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For large clusters, cellular DMFT would lead to off- 
diagonal terms in the impurity action written in the ba- 
sis of cluster momenta. The hybridization function in 
Eq. would take the form ip^ A KK , ip . However, in 
the 2 x 2 case, both in C-DMFT and DCX, the hybridiza- 
tion function is diagonal in cluster momentum. 

The DMFT mapping of the lattice model onto a pla- 
quette in a medium, allows us to make a connection be- 
tween this problem and the multiorbital Hubbard models 
which have been studied in connection with the orbitally 
selective Mott transitio n 73 i 74 i 75 . This is defined by a set 
of bands, each one characterized by a local density of a 
states, labeled by its cluster wave- vector. Notice however 
that the interaction among the orbitals, i.e. the Hubbard 

U term written in terms of ibl and ib v is more compli- 

— . — 1\ 

catcd than what has been treated in the literature and de- 
serves further investigations. The local density of states 
corresponding to the different bands can be obtained by 
setting U = and evaluating the non interacting Greens 
function Go corresponding to each cluster wave vector. 
This is plotted in figure [TJ 




' \ Cellular DMFT 
/ \ 



FIG. 1: (Color online) Tight-binding DOS for the three 
orbitals within DCA and C-DMFT. Notice that the tight- 
binding Hamiltonian within C-DMFT Eq. (Tl3)l contains off- 
diagonal elements therefore DOS does not contain full infor- 
mation about the non- interacting part of the Green's function 
Go (Go 1 = G- 1 + E) 



The formalism is easily extended to the superconduct- 
ing state introducing Nambu notation 



K 



CKT 

'-KJ. 



(9) 



Assuming singlet pairing, all the previous discussion car- 
ries through, with the cluster Greens functions and hy- 
bridization functions taking the 2x2 matrix form: 



G k (t) = -<T^ k (t)^(0)> 



^c(t) -S-Ki(-T) 

(10) 

Here is the anomalous component of the Green's 
function. Hybridization A K becomes a matrix as well 



A K (iw) 



Akt(^) 
(iw) 



a an] 



A%?(iu>) 
-A_Ki(-« 



(11) 



and the impurity problem is off diagonal in Nambu space. 

In cluster momentum basis (see Eq. (|8[)), which we em- 
ployed in this work on the 2 by 2 plaquette, DCA and 
C-DMFT share the same form of the impurity model, 
the only difference between the two schemes lies in the 
form of the self consistency conditions. This is dic- 
tated by the form of the non-interacting part of the 
Hamiltonian H and the region of momentum summation. 
In the DCA scheme the non-interacting Hamiltonian 
ifk is just the tight-binding energy et = — 2t(cosk x + 
cos k y ) — At' cos k x cos k y . In the self-consistency condi- 
tions Eqs. ^ and (O, the summation has to be per- 
formed only in the region of the patch corresponding to 
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each cluster momentum K— , i.e., 

&-E(( to+ S-* fa -2 + 0-iwM)". 

keK v v 7 7 

(12) 

The patches which correspond to different cluster mo- 
mentum K are thus completely decoupled in the self- 
consistency condition. Their coupling is only through 



the coulomb interaction. 

In the real space C-DMFT, we can still define "or- 
bitals" which correspond to cluster momenta K (sec the 
form of local quantities in Eq.©), however, these or- 
bitals are coupled through both the Coulomb repulsion 
U and the non-interacting Hamiltonian. which takes the 
following form 
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where we defined 

e£ = —t(2 + cos k x + cos k y ) — i'(l + cos cos /c y ) 
e k = ^(cos /e x — cos k y ) + t'(l + cos k x cos fc^) 
e£ = — t(cos k x — cos fcj/) + t'(l + cos fc^ cos fc y ) 

= t(2 + cos k x + cos fc y ) — i'(l + cos fc^ cos k y ) 
vl = t' sin k x sin fcj, 
4 = sin k x (t + t' cos k y ) 
vl = sin k y (t + <' cos k x ) 
vl = sin fc a (t — t' cos A: y ) 

vl = sin k y (t — t 1 cos k y ). (14) 

The unit of distance choosen here is a = 1/2 such that 
the summation over the reduced Brillouin zone in Eqs. ([6|) 
and ([7|) simply runs over k x 6 [— tt,tt] and fc y G [— 7r,7r]. 
One can readily show that this summation leads to diag- 
onal form of local quantities. 

In figure [2] we compare the local spectral function of 
the t-J model in the two cluster schemes. Notice the sim- 
ilarities of the results in particular at low energies. The 
spectral functions in both methods have a very similar 
pseudogap. Hence in spite of quantitative differences, 
which will not be investigated systematically in this pa- 
per, the qualitative physics, which is the main focus of 
this article is present in both cluster methods. Note, how- 
ever that decoupling of orbitals in DCA method leads to 
splitting of the Hubbard band into peaks which corre- 
spond to excitations of the 2x2 cluster. These finite size 
effects are strongly reduced in C-DMFT method. 

Here we comment on some quantitative differences be- 
tween the methods. The superconducting critical tem- 
perature is highest in EDCA method and reaches the 



value - 0.036t while it drops to - 0.026t in EC-DMFT. 
When the bosonic baths is switched-off, the real space 
C-DMFT maximum critical temperature in both, the t-J 
model at J = 0.3 and the Hubbard model at U = 12i is 
around ~ O.Olt . Notice that this value is close to the 
estimations in Ref. H3 for the critical temperature of the 
Hubbard model in the thcrmodynamical limit for U = 
4t. Namely, the Hubbard model at U — At within large 
cluster DCA has T c ~ 0.023*21 If we extrapolate this 
value to large U = 12t treating Tc cx J—, Tc would drop 
to ~ 0.008t which is close to the C-DMFT result. 

The existence of a finite transition temperature, the 
trends of the superconducting transition temperature 
with doping and with the strength of the superexchangc 
interaction is a robuts propery of plaquette DMFT and is 
common to all cluster schemes. It would be interesting to 
understand the convergence properties with cluster size 
within the different cluster schemes for the t-J model, as 
was done for the Hubbard model at intermediate U in 
Ref. and in the classical limit in Ref. WH 



D. Impurity Solvers 

At the heart of the cluster DMFT is the solution of 
the impurity problem Eq. (0. In this work, we used 
two different impurity solvers both based on the expan- 
sion of the impurity action with respect to hybridization 
strength. The first is the non-crossing approximation 
(NCA) which sums up all diagrams with no crossing and 
is conveniently formulated in slave particle approach^. 
The second is the recently implemented continuous time 
Monte Carlo Method^^ (CTQMC) which numerically 
samples the same type of diagrams but sums up all dia- 
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FIG. 2: (Color online) Comparison between the local spectral 
function computed in C-DMFT and in DCA with NCA used 
as impurity solver. 
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FIG. 3: (Color online) Comparison of NCA and CTQMC 
Green's functions on imaginary axis for several doping levels. 
We used real space C-DMFT. 



grams using Monte Carlo importance sampling. 

The two impurity solvers are in good agreement with 
each other on the imaginary axis, but the first method 
allows us to obtain real frequency correlation functions 
which are unavailable in the QMC approach. Both ap- 
proaches arc well suited to study the regime of interme- 
diate temperatures and dopings, close to the tip of the 
superconducting dome, separating overdoped and under- 
doped regions, which is not easily accessible with other 
techniques. 

Both impurity solvers used here require the introduc- 
tion of the cluster eigenstates obtained by the exact diag- 
onalization of the cluster, i.e., H c i uster \m) = E^ uster \m). 
To each cluster eigenstate, a pseudoparticle 
assigned, i.e., 

\m) = 4J0) (15) 

to recast the cluster part of the action to a quadratic 
form. The constraint 

Q = E a «> am = E l TO )H = 1 ( 16 ) 



which expresses the completeness of the atomic eigen- 
base, has to be imposed. 

The original problem can be exactly expressed in terms 
of pseudoparticles a m with the only non-quadratic term 
of the converted action being the hybridization between 
the cluster and the medium 



S 'ff= f 'E^Hl- + E ™ Stm - + f dr f dr' J2 al(T)a n (T)D mnn , m ,(T-T')al(T')a m ,(r'%17) 
Jo m ° T Jo Jo mnmlnl 



mnm'n' 



denoted here by 

D mim2m . imi (iuj) = ^(F Kt ) mir „ 2 A K (icj)(F K ) m 3 mi - -(S K )mim 2 Xo 1 K(* W )(^K)m 3 m 4 (18) 



K 



where 



(F K ), 
(Sq)r 



(m|^ K |n) 
(m|Sq|n) 



(m\cKf\n) 
(m|cL K; |n) 



(19) 
(20) 



Note that the effective hybridization D combines both, 
the fermionic (A K ) and bosonic bath (xo^k) the 
total effective Weiss field felt by the cluster eigenstates 
(pseudoparticles). We used lagrange multiplier A to en- 
force the constraint (fTfJ|) . 

The continuous time quantum Monte Carlo method 



samples over the diagrams generated by expanding the 
action J D[a'a] cxp(— S c i us ter — AS) with respect to ef- 
fective hybridization AS. Here AS stands for the sec- 
ond term in Eq. [171 The probability to visit each dia- 
gram is proportional to its contribution to the partition 
function which is computed by explicit evaluation of the 
cluster trace J D[a) a]e~ Scl,lst ' r (— AS) k /k\ keeping only 
single pseudoparticle in the system at each moment in 
imaginary time. In this way, the constraint Q = 1 is ex- 
plicitely taken into account. For more details, see Ref . l58l. 

In the diagrammatic method, the constraint Q = 1 is 
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imposed by letting the Lagrange multiplier A approach 
infinity. The physical observable can then be computed 
using the Abrikosov's trick^ (A)q = i = lim^oo ^f^ 1 ■ 

The coupling of the cluster to the medium, which simu- 
lates the rest of the lattice, causes the cluster eigenstates 
to decay in time. Therefore their spectral functions carry 
nontrivial frequency dependence and important informa- 
tion about various physical processes such as the RKKY 
interactions, the Kondo effect and d-wave superconduc- 
tivity. The corresponding pseudoparticlc Green's func- 
tion can be written in the form 

Gmn(u) — {id + A — Eduster — S) 1 mn (21) 

where {E c i uster ) mn = {E cluster ) m 5 nm is the energy of the 
cluster eigenstate and A is Lagrange multiplier which will 
be sent to infinity at the end of the calculation. 

Although hybridization is a small quantity compared 
to other scales in the problem, the perturbation is sin- 
gular in the sense that at zero temperature an infinite 
number of diagrams substantially contribute to the solu- 
tion of the problem. In Ref . [HI we showed a histogram (a 
distribution of the perturbation order) which is peaked 
around (E^ ra )/T where Eki n is average kinetic energy 



and T is temperature. An infinite resummation of dia- 
grams is thus necessary and the non-crossing diagrams 
are simplest to compute. 

Just like in the single site Anderson and Kondo impu- 
rity proble m 59 ' 61 , the non-crossing approximation works 
well down to some breakdown temperature, which is 
slightly below the superconducting transition tempera- 
ture. Although NCA is not exact, this approximation 
has the virtue of yielding directly real frequency infor- 
mation. In Fig. [3] we present a typical comparison of 
the two impurity solvers on imaginary axis for the clus- 
ter Green's functions of the t-J model in normal state 
close to T c . This comparison illustrates the degree of 
agreement within the two solvers on the imaginary axis. 
Notice that all the qualitative features of the evolution of 
the Greens functions with doping are seen in both meth- 
ods. Therefore we will use in this work, the strategy of 
combining information from different solvers, in order to 
draw conclusions as to the physical picture contained in 
the solution of the cluster DMFT equations of the t-J 
and Hubbard model thus avoiding the difficult problem 
of analytic continuation of imaginary time QMC data. 

In the non-crossing approximation, the pseudoparticle 
self-energies are computed from 



S TO / m (iw) = T G„/„(ie) {D nmm , n i(ie - iu>) - D m , n , nm {iu> - ie)} 



(22) 



while the physical quantities like Green's function and susceptibility arc obtained by functional derivative of the NCA 
Luttinger functional with respect to the hybridization term and are given by 



(23) 



le.mnm n 



Xcf(^) = T ^2 (Sq)m'7i'G„'7i(ie)G mm '(ze - iw)(S^ Q ) nm . (24) 



ie ,mnm' n' 



The above equations can be projected to the physical subspacc Q = 1 only on the real axis. In the limit A —> oo they 
take the form 



£ m < OT (w) = ^ J d£,f(£,)G n >n(£, + |D„mm'n'(0 + An'n'nm(-£)} 
K.nn' 



(25) 



<2kM 

XqV) 



£ (F K ) m ,„,(F K t), 

mnm'n' 

(Sq)m'n' (Sq)rin 



mnm'n' 



Here we used the following notation 
1 



D{u) 
G 



2-ki 
1 

"2ri 



[D(w + iS) - D(uj - iS)} 
[G(uj + i6)-G(Lj~i6)}. 



G n . n (£ + u)G mm , (0 + G n , n (OG* mm , (£ - w) 
I 



(26) 
(27) 



The pseudoparticle quantities (Green's functions G and 
self-energies S) exponentially vanish below a certain 
(28) threshold energy (they have X-ray singularity) which can 
be interpreted as the effective energy of the many-body 

(29) 
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state associated with the pseudoparticle. These thresh- 
olds can be removed by defining new quantities without 
threshold^, i.e., 

G(e) = G(e)//(-e) (30) 
E(e) = S(e)//(-e) (31) 



At zero temperature, the combination of the Fermi func- 
tions — jrzj^ ^ = *^f(u) ^ * s ec l ua l to unity in the 
interval [mm(0,w),mo3;(0,w)] and zero outside. 

These equations relate physical observables, like Gk 
and xq to the pseudoparticle spectral functions. The 
later represent coarse grained versions of the important 
many body excitations of the system including fermionic 
quasiparticles and bosonic collective modes. They have 
quantum numbers describing their spin, number of parti- 
cles, (which divided by the cluster size, give the density), 
and a coarse grained momentum. 



In this section we discuss cluster quantities. As dis- 
cussed in section HQ in both C-DMFT and DCA for- 
malism local quantities, such as cluster self-energies and 
cluster Green's function are diagonal in the cluster mo- 
mentum basis. Consequenctly, the physical behavior of 
the system within the cluster DMFT approach on a pla- 
quette can be summarized in the four cluster quantities 
Soo, Stto, and X^, corresponding to the eigenval- 
ues of the matrix containing onsite, nearest-neighbor and 
next-nearest neighbor cluster self-energy introduced in 
reference H^. These cluster self energies in the cluster 
momentum basis should not be interpreted as the lattice 
self energies evaluated at 4 momentum points. 

In the next few figures, we present low temperature 
self-energies for the t-J model on the imaginary axis ob- 
tained using CTQMC impurity solver. Figure 2] contains 
the data in the normal state and Fig. O same quantities 
deep in the superconductings state. 

Starting from the low temperature T = O.Oli normal 
state solution shown in Fig. 2] one notices large momcn- 



Using these quantities we can rewrite the NCA equations 
as: 



(32) 



(33) 
(34) 



Relating several experimental observables such as pho- 
tocmission spectra, tunneling spectra, and optical spec- 
tra, to the same set of pseudoparticle spectral functions, 
gives additional insights into the important excitations 
of the system. 



III. CLUSTER ONE PARTICLE GREENS 
FUNCTION, CLUSTER SELF ENERGY AND 
SCATTERING RATE 



I 

turn differentiation at small doping. The three orbitals 
evolve very differently with changes in doping and tem- 
perature as we will show in the following. The (0, 0) com- 
ponent has Fermi-liquid frequency dependence with rel- 
atively small scattering rate at zero frequency and small 
monotonic decrease of the real part of the self-energy with 
increasing doping. The two degenerate orbitals (0, 7r) and 
(it, 0) arc distinctly different from the (0, 0) orbital. The 
scattering rate around optimal doping 0.12 ~ £ ~ 0.22 
remains large (of order unity) even below the transition 
to superconducting state T ~ 0.01£. We notice in passing 
that it becomes increasingly difficult to converge the C- 
DMFT equations in the metastable normal state around 
optimal doping. Critical slowing down is observed which 
might be a signature of local (cluster) instability which 
might occur at zero temperature and might even pre- 
clude the continuation of a translationally invariant nor- 
mal state solution down to zero temperature. 

The (irir) self-energy in Fig. 2] is far the largest among 
all four and, except at very small and very large doping, 



S m ' m (cd) — y ] / d£ — — — —r Gn'n(£) | Dnmm'n' (£ — Ul) + D m > n ' nm (u! — £)| 

K,nn' J A - W J 

ImG K M = -7T ( pK W (F Kt W / J0P^ G n . n {£)G mm , (C - u>) 

mnm'n' 



mnm' n' 
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FIG. 4: (Color online) C-DMFT cluster self-energies of the t-J model using CTQMC as the impurity solver. Temperature 
T — O.Oli and system is in the normal state. Notice that Eqo is Fermi-liquid like (imaginary part vanishes at zero matsubara 
frequency below the coherence temperature) in the whole range of doping, E-n-o is Fermi liquid in the overdoped and underdoped 
regime while the scattering rate remains of the order of unity in the optimally doped regime. Finally, E^ is far the largest 
self-energy. Its real part is so large that the orbital is gapped in all doping range considered. The scattering rate is enormous 
and a pole appears on the real axis around S — 0.1. The pole is above Ef at very small doping, crosses Ef at 5 = 0.1 and goes 
below Ef for optimally doped and overdoped regime. This causes a sign change of the real part of E^. The (irir) orbital is 
thus in the Mott insulating state in most of the doping range considered. 



it does not show any signature of coherence. At 6 ~ 0.1 
it has a clear pole at zero frequency. From the above 
plot we can see that a pole is on the real axis and it is 
located above the Fermi level at small doping and crosses 
Fermi level around 5 ~ 0.1 and finally gets negative in 
the optimal and overdoped regime. A very sharp pole on 
the real axis described above is indeed confirmed by the 
NCA calculation. The consequence of the pole in self- 
energy is appearance of zeros of the Green's function as 
discussed in Ref . |43|. Physically it means that some states 
in momentum space are damped and gapped even at very 
low temperature. Figure \5\ at lower temperatures show 
that this behavior persist to temperatures much lower 
than T c . Hence even in superconducting state the large 
luttinger Fermi-surface is not recovered. The antinodal 
fermions are strongly damped and gapped even in the 
superconducting state. This is related to the occurrence 
of Fermi arcs and lines of zeros of the Greens functions* 
as noticed in Ref. 43 44.81. This phenomena was first 
noticed microscopic studies of coupled ladder a 85 i 86 and 
related proposals also appeared in recent phenomenolog- 
ical models of high Tc' o 82 ' 83 ' 84 . However, in these studies, 
the location of the lines of zeros is tied to the unklapp 
surface, while in the cluster DMFT the lines of zeros is 
a dynamical object which evolves in a highly non trivial 
way with doping. 

The pole in (7T7t) self-energy crosses the Fermi level at 
a critical doping (or at least becomes very large at low 



energies) that we denote by S^. The existence of a pole 
in the self energy appears also in Hubbard model, with 
an important difference. In the Hubbard model, the pole 
is always below Ep and therefore this "critical" doping 
51 is zero. We also want to mention that at small U = 6 
(below the Mott transition of the undoped system) in 
the Hubbard model the above mentioned pole seems to 
be absent or at least substantially reduced. This sub- 
stantiates the idea that the lines of zeros in the Green 
function appear only above a critical coupling. 

Fig. [5] shows cluster self-energies at a lower tempera- 
ture, i.e., T — 0.005i in the superconducting state. The 
(0, 0) orbital does not change very dramatically execpt 
that becomes more coherent. On the other hand, the 
(it, 0) orbital does show a dramatic effect. The huge 
scattering rate is now replaced by the large anomalous 
component of the self-energy while the scattering rate 
is severely reduced. The peak in anomalous self-energy 
seems to track T c and coincides with the point of max- 
imal scattering rate in the normal state. We will call 
this doping 5% since it corresponds to the avoided critical 
point identified in Ref. [50l. Finally, the (tt, tt) component 
of the self-energy sharpens with reducing the tempera- 
ture and the pole at ~ 0.1 is even more apparent. 
This result is quite surprising because the superconduct- 
ing state would have been expected to be more coherent. 
As we show above, coherence is only restored in three of 
the four orbitals, while the (7T7t) orbital remains gapped. 
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FIG. 5: (Color online) Same as in Fig. [4]but at lower temperature T — 0.005i in the superconducting state. The bottom row 
shows the anomalouse self-energy. The (0, 0) orbital barely changes in superconducting state. On the other hand, the large 
scattering rate in (n, 0) orbital is severly reduced in the superconducting state and the orbital becomes Fermi liquid like. Large 
scattering rate in the normal state is now replaced by a large anomalous component of self-energy (peaked around 8 ~ 0.15 
- see Fig. I18[) . Finally, the pole in (tt, 7t) self-energy sharpens and the orbital remains Mott insulating in most of the doping 
range considered. A pole of the cluster self-energy is accompanied by a line of zeros of the Green's functions in certain parts 



of the momentum space 



and persist in the superconducting state. 



Hence the Fermi surface underlying the normal state does 
not contain the Luttingcr volume at small doping. 

In Fig. [S] and [JJ we show the cluster Green's functions 
of the t-J model in the normal state and the supercon- 
ducting state at lower temperatures. The cluster Green's 
functions describe a coarse grained average of the lattice 
Green's function over some parts of the Brillouin zone. It 
is evident from Fig.[S]that the (tt, 0) orbital contains most 
of the spectral weigh (largest imaginary part of G^q) over 
the whole doping regime considered here. The (jr, tt) or- 
bital is clearly gapped since the real part of self-energy is 
too big to pick up any states inside the band as was previ- 
ously observed in the extended DMFT study of the same 
mode l 60 i 61 . The important message is contained in the 
real part of (tt, 0) cluster Green's function. The real part 
measures the particle hole asymmetry of the orbital. It 
would vanishe if the orbital is perfectly particle-hole sym- 
metric. As one can see in Figs. [Hand [7Jthe (tt, 0) orbital 
has "more weight" below Ep in underdoped regime and 
more weigh above Ep in overdoped regime. Remarkably 
it becomes almost particle-hole symmetric in the region 
of optimal doping. The exact point of particlc-holc sym- 



metry is close to ~ 0.18 which is just slightly above the 
point of maximal T c and maximal anomalous self-energy. 
Fig. [JJ demonstrates that this remarkable symmetry per- 
sist even in the superconducting state where the gap ap- 
pears in all the orbitals. 

We now compare the previous findings with the corre- 
sponding quantities in the Hubbard model displayed in 
figures [51 [11 [UJ and [JJ The Hubbard model at U = 12t 
has roughly the same superexchange as the one used 
in the previous study of the t-J model J ~ 0.3i and 
therefore we expect similar physical behaviour. We will 
demonstrate below that indeed this is the case, and high- 
light some quantitative differences between the two mod- 
els, such as the numerical values of the critical dopings 
for cluster quantities 81 and 5% ■ 

Fig. [SJ shows the four cluster Green's functions at 
T = O.Oli in the normal state. When the off-diagonal 
long range order is allowed, the system starts to develop 
anomalous components in the optimal doped regime at 
this temperature just like in the t-J model at the same 
temperature. In the metastable normal state, the (0, 0) 
orbital is again most coherent orbital and not very sen- 
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FIG. 6: (Color online) C-DMFT cluster Green's functions at T = O.Olt in the normal state of the t-J model obtained by 
CTQMC The real part of the Green's function vanishes for particle-hole symmetric situation while its positive when spectral 
weight below Ef has "largest weight" and vice versa. The (0,0) orbital does not change much with doping and remains close 
to half-filling. The (it, 0) orbital gives most of the weight at the Fermi level (has largest imaginary part at zero frequency) and 
remarkably becomes particle-hole symmetric at the doping level slightly larger than the optimally doped level (S = 0.18). The 
(tt, 7r) orbital is gapped for all doping levels. 



k=(0fi) fc=(7r,0) k=(n,n) 
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FIG. 7: (Color online) C-DMFT cluster Green's functions of the t-J model at T = 0.005t in the superconducting state obtained 
by CTQMC. Superconducting gap opens in particular in (n, 0) orbital. Particle-hole symmetry of this orbital is again evident 
from the real part of the Green's function being close to zero around optimal doping (blue curve with triangles pointing right). 
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FIG. 8: (Color online) Hubbard model cluster self-energies in C-DMFT obtained by CTQMC at T = O.Olt and U = 12t in the 
normal state. Just like in the t-J model, the (0, 0) orbital is Fermi-liquid like in the whole doping regime while (it, 0) orbital is 
coherent only in the underdoped and overdoped regime. At optimal doping (in the Hubbard model optimal doping is around 
5 ~ 0.1) the scattering rate is largest. The important difference appears in the (n,Tr) orbital. The (n, n) self-energy is the 
largest self-energy of the system just like in the t-J model. Contrary to the t-J model, the pole in the (it, it) self-energy on the 
real axis, which appeares in the t-J model around S = 0.1, is now at zero doping. The self-energy of the (7r,7r) orbital thus 
monotonically grows when approaching the Mott insulator. 
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FIG. 9: (Color online) Similar than Fig.[8]but at lower temperature T = 0.005£ in the superconducting state. Just like in the t-J 
model, the (it, 0) orbital, which is representative of the nodal part of the self-energy, becomes coherent in the superconducting 
state and the anomalous self-energy is largest around 8 ~ 0.1 where the scattering rate is largest in the normal state. The 
(tt, tt) self-energy sharpens with decreasing temperature just like in the Hubbard model showing that this orbital is in the 
Mott-insulating state in the underdoped and optimally doped regime. 



sitive to doping. On the other hand, the (tt, 0) orbital 
is clearly coherent for small and large doping and the 
scattering rate around S ~ 0.1 is of the order unity. The 
point of maximum scattering rate and maximum anoma- 
lous self-energy in the Hubbard model is however slightly 
shifted towards lower doping, (relative to the t-J model) 
i.e., 5% ~ 0.1. 

The (tt, tt) orbital is again the one with far the largest 
self-energy and scatering rate. In the Hubbard model the 
pole on the real axis crosses zero exactly at zero doping, 
hence 51 = 0. However even at optimal doping ^0.1 the 
real part of the self-energy is so large that the orbital is 
almost completely gapped. 

Fig. [9] demonstrates that the pole in E W7r does not dis- 
apear in the superconducting state. This was also the 
case in the t-J model, and it is therefore a robust feature 
of the approach to the Mott insulator within CDMFT. 
The physical interpretation is that part of the underly- 
ing Fermi-surface remains gapped even in the supercon- 
ducting state. The (tt, 0) orbital becomes coherent when 
entering the superconducting state. Its imaginary part, 
at low frequencies is maximal around S^- 

The cluster self energies in the cluster site represen- 
tation contain useful information about the range. For 



example, it has been argued that near the Mott insula- 
tor, they become long ranged, while the cluster cumulant 
remains short range d 43 ' 44 . In Fig. [T0]we show the onsitc, 
nearest neighbor and next-nearest neighbor self-energy, 
the actual output of the C-DMFT scheme. These are 
related to the eigenvalues shown above through the fol- 
lowing linear relation 
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On the heavily overdoped side of the Hubbard model 
5 > 0.16 presented in Fig. [TU] it is clear that the only 
relevant quantity is the on-site self-energy which justifies 
use of the single site DMFT in the overdoped site of the 
system. In the underdoped regime, however, the nearest 
neighbor as well as next-nearest neighbor self-energies 
are large and give rise to qualitatively different results 
than those of a single site DMFT. They rcnormalize the 



15 




iuj iuj 



FIG. 10: (Color online) The on-site and short range self-energy of the Hubbard model in the superconducting state at 
T = 0.0054. The on-site self-energy is the largest and its imaginary part vanishes for all finite dopings. The reason is that the 
pole in the (tv, tt) self-energy is now at zero doping. The non-local components of the self-energy vanish rather rapidly with 
doping. 



k=(Q,Q) ft=(7T,Q) fc=(7T,7r) 
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FIG. 11: (Color online) Hubbard model Green's function in normal state at T = O.Olt. The physics is the same as in Fig. [6] 
for the t-J model. The only difference is that the particle-hole symmetric point of the (n, 0) orbital appears around 8 ~ 0.12 
which is again slightly above the optimal doped level. 
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FIG. 12: (Color online) Similar than Fig. [TT] but for lower temperature T — 0.005£ in the superconducting state. 
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nearest neighbor and next-nearest neighbor hopping and 
induce a substantial next nearest neighbor hopping even 
for the model with vanishing bare t' . Furthermore they 
distort the Fermi surface and cause variation of coherence 
across the Fermi surface as we will show below. 

Finally, the cluster Green's functions for the Hubbard 
model are shown for two temperatures T = 0.01 and 
T = 0.005 in Figs. ED] and [Hat U = I2t. Again we no- 
tice qualitatively similar behaviour than those found in 
the t-J model. The (tt, tt) orbital is gapped in both nor- 
mal and superconducting state. The (tt, 0) orbital con- 
tains most of the spectral weight and becomes particle- 
hole symmetric slightly above optimal doping around, 
i.e., around 5 « 0.12. This particle hole-symmetry per- 
sists in the superconducting state. 

We now turn to the real frequency information. In 
most of what follows, we show results for the t-J mode, 
except when explicitly stated otherwise. 

Figure [13] shows the evolution of the CDMFT cluster 
spectral functions as a function of frequency for few dop- 
ing levels. Notice that due to symmetry (tt, 0) and (Q,tt) 
spectral functions coincide. At zero doping (not shown) 
all four orbitals are half-filled and the system is in Mott 
insulating state. 

Upon doping the system, the (tt, tt) orbital gets emp- 
tied first but in a very unusual way. Although its oc- 
cupancy gets much smaller than unity and therefore one 
would naively expect large number of hole carriers in this 
band, it remains basically gapped for arbitrary doping as 
we have established above on the basis of the CTQMC 
results. This is very unusual since one naively expects 
the orbital to be gapped only at an integer filling. Only 
at very large doping S > 0.3 the self-energy of this or- 
bital approaches the other three self-energies so that the 
self-energy becomes momentum independent and there- 
fore local. At this large doping, the (tt, tt) orbital is es- 
sentially empty and we can think of this orbital as an 
Anderson impurity model in the empty orbital regime. 

The (0, 0) orbital is also very inert in the whole doping 
range. Its density of states at the Fermi level is small 
while its occupancy only slightly decreases with increas- 
ing doping. The orbital remains close to half-filling with 
very small number of charge carriers induced in this band. 

Finally, the the (0,7r) (and (7r,0)) components have 
sharp spectral features with very strong doping depen- 
dence. In going from S = 0.3 to J = 0.1 we observe 
the narrowing of the quasiparticle width reminiscent of 
the single site DMFT, however a qualitative feature of 
CDMFT is that at smaller dopings this narrowing of the 
width is arrested, as a result of the presence of exchange 
effects as seen in slave boson studies^ and in the large 
N limit of the t-J model^. 

At low doping, the spectral function develops a pseudo- 
gap on the scale of J with most of the coherent spectral 
weight below the Fermi level and a small fraction of it 
above the Fermi level. This is a general feature of the 
approach to the Mott transition in cluster DMFT and 
has been seen in earlier studie a 41 ' 42 i 45 i 60 i 61 . 
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FIG. 13: (Color online) Evolution of the EC-DMFT cluster 
spectral functions of the t-J model with doping in normal 
state at T — 0.025t ~ T c using NCA as an impurity solver. 
The upper panel shows the spectra in the interval between 
[—St, 4t] where Hubbard band is clearly observed. The lower 
panel shows the region near the Fermi level. 




FIG. 14: (Color online) EC-DMFT cluster self energies on 
real axis in normal state of the t-J model computed by NCA 
at T = 0.025rj ~ T c . 



The important message contained in figure [T3]is that 
the momentum differentiation at small doping is very 
large. The (tttt) orbital remains gapped at all dopings. 
It is in the Mott insulating state at low doping and be- 
comes empty in the overdoped regime, hence it under- 
goes a band insulator to Mott insulator transition with 
decreasing doping. Most of the dynamical information of 
the active degrees of freedom representing the electrons 
close to the Fermi surface of the lattice model is however 
contained in the (0, tt) and (tt, 0) components. 

The frequency dependence of the cluster (tt, 0) self- 
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energy and its evolution with doping on real axis is shown 
in Fig. [T31 At small doping, the hole-like scattering rate 
{ijj < 0) is large while the electron-like (uj > 0) is small. 
Around optimal doping, the self-energy is roughly linear 
in frequency, however, with large zero-frequency value. 
In this regime, there is still a large particle-hole asym- 
metry in the scattering rate. While the hole-like part 
is linear with relatively small slope down to lu = —0.5t, 
the electron part is increasing only in the small region 
up to lu = 0.15t with larger slope. Only in the strongly 
overdopcd system, the self-energy becomes roughly par- 
ticle hole symmetric at low frequency. This particle-hole 
asymmetry in scattering rate can be contrasted with the 
approximate particle-hole symmetry in the one-particle 
green's function at optimal doping. The combination of 
the real-part of the self-energy and the band-structure, 
leads to approximate recovery of this symmetry in the 
local one-particle spectra at optimal doping although it 
is absent in the scattering rate. 

A. Zero frequency Quantities 

Further insight into the nature of the CDMFT solution 
can be obtained by examining the cluster self-energies at 
zero frequency. In Fig. Ql}] we display the CTQMC self- 
energy for the t-J model at the lowest Matsubara fre- 
quency as a function of doping. In the overdoped side, 
the real parts of all four self-energies merge therefore the 
self-energy becomes local. The single site DMFT is ade- 
quate. The coherent quasiparticle peak at the Fermi level 
is formed and arrises mainly from the (vr, 0) orbital. The 
reason is that the non-interacting density of states for 
(0,0), (tt, tt) and (0, tt) orbital extends roughly between 
[-44,0], [0,4i] and [-2t,2t], respectively (see Fig. [J). For 
the momentum independent self-energy, the Friedel sum- 
rule dictates that the effective chemical potential \i— E(0) 
is at the corresponding non-interacting chemical poten- 
tial fiQ which is slightly below zero frequency. The (tt, tt) 
orbital is therefore empty being a band insulator-like. 
At smaller doping, this orbital acquires enormous real 
part of E'(0) which pushes effective chemical potential 
/i — E'(0) far below the band edge of the tight-binding 
Hamiltonian. This orbital is therefore in the Mott in- 
sulating state for smaller dopings. The insulating state 
in this orbital does changes the nature from band-like to 
Mott like insulator. 

For the (0, 0) orbital, the effective chemical potential is 
close to its upper band edge. The non-interacting density 
of states at the band edge for this orbital is small (see 
Fig. [T]) and only a very small number of charge carriers 
are doped into the orbital. Therefore it remains close to 
Mott insulating state with small scattering rate at the 
Fermi level. 

The {tt, 0) orbital is slightly less than half-filled in the 
doping range considered here, and the real part of self- 
energy smoothly increases with doping (see Fig. [15)) such 
that the effective chemical potential \i — E'(0) is positive 
in the underdoped side (carriers are hole like) and neg- 
ative in the overdoped side (carriers are electron like). 
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FIG. 15: (Color online) C-DMFT cluster self energies at the 
lowest Matsubara point E(iojo) ~ E(w = 0) versus doping in 
CTQMC at T = O.Olt in normal state of the t-J model. 
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FIG. 16: (Color online) Phase shift (see definition in the text) 
versus doping for the four orbitals of the cluster obtained by 
CTQMC within C-DMFT method in the t-J model. 



Close to the optimal doping, effective chemical poten- 
tial is close to zero which makes orbital approximately 
particle- hole symmetric at low frequency (see Fig[T51). 

The (tt, tt) self-energy acquires a pole on real axis 
around Sf ~ 10% doping which can be identified in FigfTSl 
as a divergent point of E ffT (0) and zero of E^.„.(0). 

Figure [16] describes the low energy phase shift in each 
orbital. It is defined by Sk = arg(GK(i0+)). Phase 
shifts are defined modulo tt. Notice two important fea- 
tures: at very small doping the phase shifts in all channels 
are close to zero. They confirm the picture suggested 
in Ref. in which the cluster degrees of freedom are 
weakly affected by the surroundings. It is reminiscent 
of the RKKY phase of the two impurity Kondo model. 
The system reaches the unitarity limit, as the phase shift 
crosses tt/2 in the (0,7r) channel near 6 = 0.18. This is 
because the real part of the cluster Green's function at 
zero frequency vanishes around optimal doping as shown 
in Figs. [5] and [7] 

The indication for the existence of an anomaly around 
62 c is seen most clearly in the imaginary part of the real 
frequency electron self-energy at (0, tt) evaluated at zero 
frequency. We display EDCA-NCA calculations of it in 
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FIG. 17: (Color online) upper panel: The imaginary part of 
the cluster (n, 0) self-energy at zero frequency as a function 
of temperature for few doping levels. Computed with EDCA 
using NCA as the impurity solver. 

lower panel: An estimation of the coherent scale in the normal 
state of the t- J model (black dots) and transition temperature 
to superconducting state (red dots) within EDCA using NCA 
as the impurity solver. 



the upper panel of Fig. [T7] . At large and small doping 
the scattering rate is small as expected for a Fermi liquid. 
Remarkably it becomes very large in the region near op- 
timal doping when the critical temperature is maximal. 
This doping point corresponds to the above defined crit- 
ical doping 62 which is in NCA around 0.18. The tran- 
sition to the superconducting state severely reduces the 
scattering rate eliminating the traces of the underlying 
critical behavior. A coherence scale, estimated from the 
scattering rate, is plotted in the bottom panel of Fig. [17] 
and shows it tends to vanish close to the point of maximal 
superconducting transition temperature. 

Figure fTS] show the CTQMC results for the scatter- 
ing rate within CDMFT and confirms the incoherence 
of the optimally doped system. The imaginary part of 
the self-energy at the first Matsubara point is small for 
both, the undcrdoped and overdoped system, while it is 
peaked at optimal doping. The peak is slightly shifted 
with temperature and, if the normal state is continued 
below the superconducting transition temperature, the 
peak of scattering rate coincides with the maximum of 
the anomalous self-energy which traces maximum of the 
transition temperature (see Fig. IT^j) . The scattering rate 
is severely reduced in the superconducting state when off- 
diagonal long-range order is allowed in the calculation. 
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FIG. 18: (Color online) upper panel: Imaginary part of (7r, 0) 
cluster self-energy at the lowest Matsubara frequency inT 
versus doping for three different temperatures obtained by 
CTQMC for the t-J model. The scattering rate is peaked at 
optimal doping. 

lower panel: The large imaginary part of self-energy of the 
normal state (black curve with circles) is severely reduced in 
superconducting state (green curve with squares). The scat- 
tering rate is peaked at the point of maximal anomalous self- 
energy (blue curve with diamonds) tracking the point of the 
highest Tc- 



IV. SUPERCONDUCTIVITY, TUNNELING 
DENSITY OF STATES, FERMI ARCS AND 
NODAL QUASIPARTICLES 

The superconducting state is characterized by an order 
parameter (ck-fC-kj.) = J~v.{t = 0) and by the presence of 
a frequency dependent anomalous component to the self 
energy. In Nambu notation, the self-energy in (jt, 0) and 
(0, 7r) orbital takes the following form 



E k (mj) = 

and the corresponding Green's functions is 



(39) 



<2k(«w) 



(40) 



The sign of the anomalous components chosen by the 
system is E£g = -Eg™. Within C-DMFT, this is pre- 
cisely the nearest neighbor self-energy and its lattice ana- 
log (using the original C-DMFT periodizationiS) takes 
the form Ek = ^(cosfc^ — cosfc y )Eg™. 

The anomalous self-energy E^q is plotted in Fig. [TIB 
The upper part of the figure shows the CTQMC re- 
sults within C-DMFT while the lower part shows the 
NCA-results within EDCA. In both cases, the function is 
monotonically decreasing with imaginary frequency and 
is largest at optimal doping. Furthermore, at the low val- 
ues of the Matsubara frequency the anomalous self energy 
exhibits a fast upturn and sublinear frequency behavior 
that gets less pronounced as the doping is reduced. This 
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FIG. 19: (Color online) upper panel: C-DMFT t-J model 
anomalous self-energy computed by CTQMC at T = 0.0054 ~ 
0.5T™ ax . 

lower panel: t-J model anomalous self-energy within EDCA 
computed by NCA at T = 0.033* ~ 0.9T™ ax 



trend is likely due to the reduction of density of states in 
the pseudogap region. 

The anomalous self energy obeys a spectral represen- 
tation 



dwImS™(w) 



7T IU)„ 



The infinite frequency value of the self energy van- 
ishes in the Hubbard model but is non zero in the t — J 
model and is related to the order parameter of the system 
JF q (r = 0) through the following exact relation: 



sr(oo) = - 



(i 



S) 2 ^ 



Jk- q ^ q (r = 0). (41) 



where _F q (r = 0) = (c qT c_ qi ). 

Notice that simpler mean field theories of the t-J model 
such as the slave boson mean field theory^ assume only 
the static, frequency independent anomalous self energy. 
Other approaches based on the equation of motion for 
the Hubbard operators^ capture a frequency dependent 
order parameter but neglect the static infinite frequency 
component. A similar analysis of the pairing interaction, 
has recenty been carried out for the ladders in Ref. [9l|. 

The existence of a finite value of the anomalous self en- 
ergy of the t-J model at infinite frequency, should be in- 
terpreted as the existence of a non zero value for anoma- 
lous self energy in the Hubbard model at a scale of the 
order U. 




FIG. 20: (Color online) Order parameter in C-DMFT com- 
puted with CTQMC at T = 0.5T cmax . The critical tempera- 
ture (in units of t) for few doping values is also displayed. 



The value of the anomalous self-energy at zero- 
frequency and low temperature and the gap (defined as 
the distance between the positive and negative energy 
peaks in the tunnelling density of states divided by two) 
is similar in all versions of the cluster DMFT. For the 
parameters used in our study, (J/t = 0.3, near optimal 
doping) the anomalous self energy is of the order of unity 
at low temperature (see the upper panel of Fig. [T9|) . 

On the other hand, Tc , the superconducting order pa- 
rameter, and the value of anomalous S an (oo) are more 
sensitive quantities, and differ between the various clus- 
ter schemes. The schemes with higher Tc (extended ver- 
sions of CDMFT) show slower decrease of the anoma- 
lous self-energy, larger infinite frequency component of 
the anomalous self energy and larger value of the super- 
conducting order parameter. In C-DMFT the maximum 
value of the order parameter is around 0.02 (sec Fig. [2TJ)l 
which is approximately 8 times smaller than maximum 
achieved in EDCA . Consequently, the static pairing in 
C-DMFT is very small while it reaches almost 1/3 in 
extended versions of the cluster DMFT (both in EDCA 
and and EC-DMFT), i.e., the magnitude of the anoma- 
lous self-energy at infinity as compared to the value at 
zero shown in Fig. 1191 

From the anomalous Greens function we can extract 
the order parameter, i.e., the anomalous Green's func- 
tion at equal time T^q(t = 0). The order parameter 
versus doping as obtained by the CTQMC and C-DMFT 
is shown if Fig. [20j It has a dome-like shape and tracks 
the value of the critical temperature, just like in BCS 
theory. In the same figure, we also display critical tem- 
perature Tc at optimal doping. Due to critical slowing 
down in the region of transition, many DMFT iterations 
arc needed to determine the critical temperature. 

The temperature dependence of the related quantity, 
the anomalous self-energy at infinity, computed with 
NCA is shown in Fig. It has a clear BCS-like tem- 
perature dependence saturating to a value of the order 
of rsj 0.3 which is around 1/3 of the zero frequency value. 

Using NCA we can we can examine directly X (uj) and 
E (u>) on real axis. There are several features in the fre- 
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FIG. 21: (Color online) Anomalous self energy T, an (tu = oo) 
as a function of temperature within EDCA using NCA as the 
impurity solver at optimal doping 5 — 0.18. 

quency dependence of the anomalous self energies de- 
picted in figure f2U which exhibit noticeable departures 
from the standard Migdal Eliashberg theory of super- 
conductivity. First, the real part of the self energy does 
not change sign. This indicates that the interaction is 
attractive over the whole frequency range. There is no 
characteristic energy corresponding to WDebye where the 
interaction turns from attractive to repulsive. Further- 
more, the spectral function displays significant spectral 
weight not only at a scale of order J, but also at the scale 
of order t, extending all the way to very high frequencies. 
Several scales can be clearly identified in the anomalous 
self-energy: the size of the SC gap in one-particle spectra 
~ O.lt (see Fig. [25] and the discussion of the figure later 
in this section ), the spin exchange J, the hopping t and 
a scale of the order of half the bandwidth ~ 3t. 

It is useful to momentum resolve the one particle spec- 
tra at low energies to understand the origin of the low 
energy quasiparticle excitations in the system. This re- 
quires the choice of a periodization scheme. For simplic- 
ity we use the cumulant periodization scheme introduced 
in reference I2lll43ll44l . A more detailed discussion of the 
periodization problem will be given elsewhere^. Here we 
focus on the temperature dependence, which require the 
finite temperature techniques described in this paper. 

The results are shown in Fig. [23] As shown in earlier 
wor k 47 ' 48 , C-DMFT is able to produce Fermi arcs in the 
nodal region. The advantage of the CTQMC technique 
relative to other solvers, is that it allows to investigate, 
for the first time within CDMFT, the temperature de- 
pendence of the arcs. As shown in Fig. [22] the Fermi arcs 
shrink with decreasing temperature, reminiscent of recent 
experiments^ on cuprates. The physical mechanism for 
the formation of the arcs, and their shrinking with de- 
creasing temperature is the shift in the real part of the 
momentum dependent self energy, which is enhanced in 
the cumulant periodization. While the validity of this 
periodization down to zero temperature, with the conse- 
quent formation of lines of zeros and Fermi pockets, is at 
this point a conjecture that deserves further study, there 
is no question that the formation of the arc and their 
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FIG. 22: (Color online) Anomalous self energy on real axis 
within EDCA using NCA as the impurity solver. The top 
panel shows the extended region in frequency space while the 
lowest panel zooms in the low frequency part. 



temperature dependence, at finite temperatures is a ro- 
bust property of the cellular DMFT treatment and is vis- 
ible in other periodizations. Therefore, the results of this 
paper together with the earlier zero temperature results 
of Rcf. [43] are consistent, at the qualitative level, with 
both the recent De Haas Van Alven measurements^ and 
photoemission measurments^. With decreasing temper- 
ature the Fermi arcs, evolve into a small pocket at a finite 
distance from a line of zeros which darkens one side of 
the pocket. 

The arcs are increased with doping and they develop 
into a banana shape structure. The Fermi surface at 
optimal doping in superconducting state is displayed in 
Fig. [2H Notice the sharp quasiparticles in nodal region 
and gap in the antinodal region. 

We now turn to another observable, the superconduct- 
ing tunneling density of states and its doping dependence 
displayed in Fig |25l This quantitity has been extensively 
investigated experimentally^. On a broad energy scale, 
there is considerable particle hole asymmetry in those 
curves and the positive frequency part decreases as we 
undcrdoped. This is expected on very general grounds 
for a doped Mott insulator-S^^^^. 

Remarkably, around the optimal doping S ~ 0.18 
the curves are more particle hole symmetric at low fre- 
quencies. The asymmetry in the superconducting state 
evolves from the asymmetry of the underlying normal 
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FIG. 23: (Color online) Momentum resolved spectral function 
in underdoped regime (5 = 0.09) at zero frequency above and 
below Tc in (metastable) normal state. We use CTQMC and 
C-DMFT. 




FIG. 24: (Color online) Momentum resolved spectral func- 
tion in optimally doped regime in superconducting state 
(T = T c /2) at zero frequency. We use CTQMC and C-DMFT. 



state. To confirm this we plot the density of states of the 
underlying normal state with dashed lines in the lower 
panel of figure [25l It is clear from figure [25] that the 
same magnitude of the assymetry is present in the nor- 
mal state local density of states as in the superconducting 
tunneling density of states. 




FIG. 25: (Color online) Tunneling density of states (local 
spectral function) in EDCA obtained with NCA at constant 
T — 0.025t. The lower pannel is a blow up of the low en- 
ergy regions. The dotted lines show the normal state spectral 
function for the same doping and temperature. 



The low energy slope of the tunnelling density of states 
is only weakly doping dependent, as was shown in the 
CDMFT exact diagonalization study of the Hubbard 
model**. 

Besides the considerable particle hole asymmetry at 
low doping , there are several features in figure [35] which 
are in qualitative agreement with experiments 1 ^ 1 -. For 
example the dip-hump feature in the tunneling density 
of states, in the unoccupied part of the spectra. 

Another surprising aspect of the tunneling is that the 
increase in the gap with decreasing doping is correlated 
with a decrease in the intensity of the coherence peaks. 
This is the oppositve of what is expected for a BCS su- 
perconductor where the growth in coherence peaks corre- 
lates with an increase in the superconducting gap. This 
observation can also be understood in terms of the two 
gap picture. The gap in the tunnelling density of states, 
(maximum between the coherence peaks) is controlled by 
the gap originating from the normal component of the self 
energy. This gap increases with decreasing doping. On 
the other hand the degree of coherence is controlled by 
the anomalous self energy which decreases with decreas- 
ing doping as shown in fig 1251 

The two gap picture of the cuprates, has re- 
cently emerged from the analysis of numer- 
ous experiment s 102 ' 103 ' 104 ' 105 . It is also part 
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energy periodizationiii, i.e 
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FIG. 26: (Color online) upper panel: doping versus the chem- 
ical potential for the Hubbard model at T = 0.005i. It shows 
linear dependence and downturn at small doping. This could 
point to divergence of the compressibility zero doping, 
middle panel: the nodal quasiparticle residue Z versus dop- 
ing. It is slowly increasing in the underdoped and optimally 
doped system and increases rather rapidly in the overdoped 
system. 

lower parnnel: fermi momentum along the nodal direction 
versus doping. 



of various phcnomenological pictures of cupratc 
supcrconductor a 82 ' 83 i 106 . This picuture has been 
recently put on a microscopic basis by Variational 
Cluster Approach studies^ and C-DMFT studies of the 
Hubbard^-. 

In the optimally doped regime the gap value is of 
the order of A ~ 0.09t. This value was obtained from 
fig. but similar values result from analytic continua- 
tion of CTQMC data to real axis. As discussed above, 
the critical temperature Tq strongly depends on clus- 
ter scheme employed. Using the maximum Tq of each 
scheme, we can determin the ratio 2A/Tc at optimal 
doping. In EDCA Tq DCA ~ 0.036i and 2A/T C ~ 5, 
in EC-DMFT T% C - DMFT ~ 0.026< and 2A/T C ~ 7 
and Tg- DMFT ~ 0.01* therefore 2A/T C ~ 18. In con- 
ventional superconductors described by BCS theory, this 
ratio is universaly equal to 2A/Tc = 3.5 but increases 
in the strong coupling Eliashberg theory. The cluster 
DMFT superconductivity is thus in the very strong cou- 
pling limit when compared to conventional superconduc- 
tors. Recent experiments on Bi2213i2i seems to suggest 
that the ratio 2A/T C is close to 8.0 being somewhere be- 
tween the two limits of extended and non-extended ver- 
sion of the CDMFT schemes. 

In Figs. [27] and [26] we present some insights into the 
nodal quasiparticlcs of the Hubbard model as obtained 
from the CTQMC results shown in Figs. [Hand [H The 
self-energy in the nodal region is obtained from the self- 



S ( k ) = TpR=(0,0) + Sfl=(i, ) COSk x 

+ £fl=( ,i) cos k y + cos(fcx) cos(ky)T, R={1A j}i2) 

This allows us to determine the position of the Fermi 
momentum /i — e^ F — Sk F (u> = 0) = and quasiparticle 
renormalization amplitude Z = 1/(1 — dS(kir)/dw). In 
Fig-HUwe plot Z no dai and along the nodal direction in 
the superconducting state where the coherence is estab- 
lished and quasiparticles are well formed. Fermi surface 
is close to (7r/2,7r/2). The renormalization amplitude Z 
is very slowly increasing with doping in the underdoped 
and optimally doped regime but has a fast upturn once 
the normal state gets more Fermi-liquid like. 

The evolution of the nodal velocities at very low tem- 
peratures and its consequences for the superconducting 
gap in the nodal as well as antinodal region was recently 
studies in Ref. [Hi] using exact diagonalization as the im- 
purity solver. Here we confirm all the qualitative trends 
in the doping dependence of these quantities, using the 
CTQMC solver. 

In Fig. [23 we plot quasiparticle Fermi velocity per- 
pendicular to the Fermi surface and anomalous velocity 
parallel to the Fermi surface in the nodal region. The 
velocities are defined by 



Vnodal ^nodal I "T" 



\ dk± dk j 



I'A — Z noc ial 



X^anomal 



dk 



(43) 
(44) 



It is clear from the Fig. [27] that the nodal velocity is 
almost constant in the underdoped, optimally doped 
and lightly ov erdo ped regime, compatible with obser- 
vation in Ref. 1 1081 . The anomalous velocity, however, 
is of dome-like shape and tracks the critical tempera- 
ture. The anomalous velocity measures the slope of the 
superconducting gap at the node and its downturn in 
the underdoped regime suggest that the superconduct- 
ing gap at the node decreases with decreasing doping. 
This surprising result is in accordance with recent Ra- 
man experiments^^ and angle resolved photoemission 
measurements^^ showing that the superconducting gap 
at the node in the deeply underdoped regime indeed de- 
creases. 



V. OPTICAL CONDUCTIVITY 

We now turn to the optical conductivity, which we dis- 
play in figure[55]for the t-J model. This quantity has been 
investigated both theoretically and experimentally over 
the last twenty years. For reviews see Refs. HHJlfl. The 
integrated spectral weight is a measure of the number 
of carriers, and its evolution with doping has attracted 
considerable attentioni^29^ii^±^±i^±£. 
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The starting point of a theory of the optical conduc- 
tivity is the Kubo formulae, 



FIG. 27: (Color online) upper panel: the nodal velocity per- 
pendicular to the Fermi surface versus doping. It is almost 
flat in the underdoped and optimally doped system and in- 
creases rapidly in the overdoped system. 

lower pannel: the nodal anomalous velocity parallel to the 
Fermi surface versus doping has a dome-like shape with the 
peak around optimal doping. 
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where the plasma frequency is evaluated from the expectation value of the projected kinetic energy and the operators 
c and ct are projected fermions of the t-J model. 

In principle the evaluation of the optical conductivity within CDMFT requires the evaluation of the vertex function, 
since current vertex corrections are non vanishing in plaquettc C-DMFT. However for DCA in a plaquette, we have 
shown that these corrections vanish^. This suggest, that as a first step in investigating optical conductivity we can 
neglect current vertex corrections, and evaluate the conductivity from a convolution of the Greens functions 
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Within C-DMFT, in the regime where the NCA solvers 
can be used, the f-sum rule 

roc 2 

/ a '(x)dx = ^f- J2 [d 2 £k/dk 2 a ]n k(I (47) 

^ k,o",a — (x,y ) 

is obeyed within a few percent, suggesting, that even for 
C-DMFT, where the vertex corrections are non vanish- 
ing, the corrections introduced by this effect are small. 
Notice that the right-hand side of the Eq. (|47|) is propor- 
tional to the kinetic energy of the low energy model if 
this model contains nearest neighbor hopping only. 

Formula (|46[) depends on the momentum dependent 
Greens function and therefore on the periodization 
scheme used and the cluster method employed. The qual- 
itative features discussed in this paper and the behavior 
of the integrated quantities are common to all methods. 



The optical conductivity has been modeled as cither 
a one component or a two component system via an 
extended Drude analysi o 117 ' 118 . The two component 
parametrization consists of a Drude peak and a mid in- 
frared feature. 

The cluster DMFT results for the optical conductivity 
of the t-J model are shown in figurerjHl We show the evo- 
lution of the optical conductivity with doping at various 
temperatures. In the very underdoped regime, there are 
clearly two components to the optical conductivity with 
an optical pseudogap, opening as a function of tempera- 
ture. On the other hand beyond i5 = 0.1 one can describe 
the optics in terms of one broad feature which narrows 
as the temperature is reduced. 

It is customary to parameterize the optical conductiv- 
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FIG. 28: (Color online) Optical Conductivity at different dop- 
ings and temperatures within EDCA using NCA for the t-J 
model. The optical conductivity evolves from a sum of two 
distinct features at low dopings to a broad Drude feature at 
high dopings. 
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FIG. 29: (Color online) Integrated spectral weight in normal 
state as a function of temperature in underdoped regime for 
few different cutoff frequencies. EDCA and NCA was used. 



ity in terms of a generalized Drude model 
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This parameterization relates the real and imaginary part 
of the optical conductivity in a given energy range to two 
functions, / \ •> and 1 
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The quantity to 2 is determined from a requirement in- 
volving the energy range in which the parameterization 
is used, namely 



u> 2 r A 



(52) 




FIG. 30: (Color online) effective mass and plasma frequency 
as a function of doping. Obtained from optical conductivity 
of Fig. l28l 



where A is the high energy cutoff. 

Figure [30] describes the evolution of the plasma fre- 
quency and effective mass versus doping in the t-J model. 
The plasma frequency vanishes at half-filling and lin- 
early increases at low doping. The optical mass is weakly 
doping dependent, and changes from approximately 3 in 
overdoped regime to 5 in underdoped regime with largest 
slope at optimal doping. Weak doping dependence of the 
effec tive mass of the same magnitude was pointed out in 
Ref. [Til 

Given a parametrization of the optical conductivity as 
a sum of a few poles, the optical mass measures the ratio 
of the total spectral weight compared to the weight in 
the zero energy pole, representing the Drude peak. If the 
transitions between the upper and lower Hubbard band of 
the Hubbard model are included in u> p , i.e., A > U, than 
Lu p is finite approaching the Mott transition and conse- 
quently the optical mass diverges. On the other hand, ex- 
cluding transitions between the Hubbard bands, results 
in lu p vanishing as the Mott transition is approached. 
In the t-J model, the upper Hubbard band is projected 
out, therefore the optical mass is always finite. As long 
as the transitions into the upper Hubbard band are ex- 
cluded, the plasma frequency u> p of cluster DMFT and 
single site DMFT are not too different. Notice, how- 
ever, that m* jm is enhanced in cluster DMFT relative 
to single site DMFT (not shown) because superexchange 
transfers optical weight from the low energy to the inter- 
mcditate energy range <~ J. 

The optical spectral weight us 2 is in general function of 
temperature and cutoff A, i.e., 



W(A,T) 



(53) 



In experiment, the cutoff is usually chosen such that the 
intcrband transitions are absent (A ~ lcV). The inter- 
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FIG. 31: (Color online) The optical conductivity a(u) of the t- 
J model is proportional to w" 2 / 3 in the intermediate frequency 
region for optimally doped system. 
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FIG. 32: (Color online) Superconducting and normal inte- 
grated spectral weight (N e ff) as a function of cutoff fre- 
quency. Optical spectral weight which collapses to delta func- 
tion in superconducting state comes from a very extended 
energy interval (~ 3t). We employed EDCA and NCA. 

band transitions or transitions into the upper Hubbard 
band are absent in the t-J model therefore this require- 
ment is taken into account automatically. 

The optical pseudogap which separates the two com- 
ponents of spectra and is seen as a dip at the scale of J in 
Fig- 1221 is quite large in the underdoped system S ~ 0.05. 
One could expect that the integral spectral weight W(A) 
for small enough A ~ J might start to decrease below a 
certain characteristic temperature of a pseudogap. How- 
ever, as shown in Fig. [2U there is no sign of such a de- 
crease for any cutoff frequency A or any temperature. 
Although the pseudogap gap clearly increases with tem- 
perature, the Drude peak more than compensates for this 
spectral weigh loss and W incresases as T decresases. 

Near optimal doping, the optical conductivity displays 
remarkable power laws in an intermediate asymptotic 
regime. These power la ws w ere first pointed out by N. 
Bontemps group in Ref. Ill9l . The power laws, and the 
possibility to a connection to an underlying quantum crit- 
icality, has been a subject of several recent experimental 
papers^. CDMFT provides a natural explanation for 
these anomalies^. These powerlaws were seen in ex- 
act diagonalization of much larger systems^, indicat- 
ing again the power of the cluster DMFT when it can 
be compared with available exact results. The power of 
the optical conductivity is very close to 2/3 as seen in 
figure [3~T1 but an analytic derivation of this result is not 
available. 
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FIG. 33: (Color online) Difference between the low tem- 
perature and 300K optical spectral weight integrated up to 
A = 6f. The cluster data are computed within EDCA for 
the t-J m odel. The single site DMFT data are reproduced 
from Ref. 1 1231 and correspond to the Hubbard model with 
the same superexchange interaction. The error bars are due 
to extrapolation of spectral weight to zero temperature from 
finite temperature results {T min ~ 0.5T C ). 



A surprising aspect of the physics of strongly correlated 
materials, is that low energy phenomena affects the spec- 
tra of the material over a very large energy scale. This 
general phenomena is illustrated in Fig. 1321 which shows 
the integral of optical spectral weight W(A) in the nor- 
mal and the superconducting state. Low energy phenom- 
ena like the onset of superconductivity which involves a 
scale of a fraction of J, involves redistribution of opti- 
cal weight of the order of At « leV which is 40 times 
more than the gap value. A theoretical insight from our 
calculation is that the high frequency redistribution of 
weight comes from the anomalous Greens function T *T 
in Eq. l|46p and hence can not be observed in the density 
of states or ARPES measurements. The large range of 
redistribution of spectral weight has been als o measured 
on cuprates and pointed out in Ref. 1 1 9ll 1 1 2t 

It is useful to compare the results for the temperature 
dependence of the integrated spectral weight of cluster 
DMFT with those of single site DMFT as reported by 
Toschi et.alii2£ . These are displayed in figure l3lfl There 
are two important observations, first the sign and the 
order of magnitude of the effect is similar in single site 
and in cluster DMFT. This indicates that local quan- 
tities can be reliably computed in this framework, and 
do not change dramatically as the sizes of the cluster is 
increased. Second, the doping dependence of this low en- 
ergy kinetic energy difference W(Q) — W(300K) has oppo- 
site slope in cluster DMFT than in single site DMFT. An 
interesting question is whether the existent experimental 
data agrees better with the singl e site or clus ter DMFT. 
It seems that the results in Ref. Il2llll24lll25l are in bet- 
ter agreement with the plaquette DMFT, but a more 
detailed comparison between theory and experiment, is 
left for future studies using the more realistic band struc- 
ture of each compound and a more precise periodization 
scheme. 
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We can also compute the temperature dependent resis- 
tivity, using EDCA in which the vertex correction vanish 




Notice that while the scattering rates at zero frequency 
tend to saturate at high temperature (see Fig. ITT)) . the 
resistivities do not as seen in Fig [Ml 

Notice that the scaling of the resistivity with the num- 
ber of holes is approximately obeyed, and that the maxi- 
mum amount of linearity is obtained near optimal doping. 
More detailed comparison with experiments will require 
a more realistic modelling of the band structure and a 
detailed investigation of the dependence of this quantity 
on periodization sheme used. 

Finally, since we have access to both the real and imag- 
inary part of the optical conductivity we can compute 
the superconducting stiffness, defined as the strength of 
the delta function peak in the superconducting state. Its 
temperature and doping dependence close to T c is dis- 
played in figure [35l In optimal and ovcrdoped regime, 
stiffness is linear function of temperature close to the 
transition while it is substantially reduced in the under- 
doped regime due to opening of the pseudogap. Si mila r 
trend was found in cuprates as pointed out in Ref. 1 1261 . 

With NCA, we are not able to reach sufficiently low 
temperatures to address the crucial issue of the doping 
dependence of the linear term of the supcrfluid stiffness. 
Instead we use the techniques of section IIVI to evaluate 
the low temperature behavior of this quantiy 

„.(0) - p s (T) = —T = e ^-^f(k B T)^l (55) 

where Z n , va , and v f were evaluated in section IIVI and 
vq is the band fermi velocity. The coefficient b can be 
evaluated directly from imaginary axis data of CTQMC 
and is plotted in Fig. [351 

CDMF captures the weak dependence of b on dop- 
ing, which was a subject of intensive experimental 
investigation s 127 ' 128 . More detailed studies of this quan- 
tity in C-DMFT, including vertex corrections, and more 




FIG. 36: (Color online) Coefficient of the linear term of the 
superconducting stiffness 7r(p s (0) — p s (T))/T versus doping 
for the Hubbard model using CTMQC. 

investigations of the periodization dependence of this 
quantity, as well as the related Big slope of the Raman 
scattcringiSi is certainly warranted. 



VI. MAGNETISM SUPERCONDUCTIVITY 
AND SPIN RESONANCE 

In this section we turn to the magnetic properties, 
starting from the cluster magnetic quantities 8 - ' 129 ' 130 . As 
in the rest of the paper, we confine ourselves to the study 
of minimal models, in this section the t-J model with 
t' = 0. Notice however, that it is known from numerous 
studies that the presence of a next nearest neighbor hop- 
ping t' affects significantly the region of stability, of the 
magnetism, and can suppress it altogether—. 

The static cluster susceptibilities are displayed in fig- 
ure 1371 These are coarse grained versions of the momen- 
tum dependent magnetic susceptibility, over the different 
regions of the Brillouin zone of the size of one quarter 
of the first Brillouin zone. While the cluster suscepti- 
bilities are relatively smooth functions of doping, they 



28 




0.2 

T/t 



FIG. 37: (Color online) The cluster spin susceptibilities of 
the t-J model versus temperature at two different dopings 
obtained by EDCA and NCA. 
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FIG. 38: (Color online) The Q = cluster spin susceptibil- 
ity versus temperature at different dopings for the t-j model 
within EDCA. 

clearly demonstrate that the spin fluctuations in differ- 
ent regions of the Brillouin zone have dramatically differ- 
ent behavior with doping and temperature. The (n, tt) 
component, dominated by the Xtttt susceptibility strongly 
increases as we approach half filling. In contrast the uni- 
form component, xoo decrease as doping is reduced, a 
signal of the opening of the pseudogap. The same is true 
of XttO- Hence, an interesting property of the pseudogap 
state is the increase of staggered magnetic fluctuations 
with the opening of the pseudogap. A similar contrast 
between the staggered and uniform response, is seen in 
their temperature dependence. We see that while the 
uniform response decreases with temperature in the un- 
derdoped regime, the staggered response increases. 

We now proceed to uniform spin susceptibility shown 
in Fig. The q = susceptibility at zero doping 

displays the characteristic behavior of the Heisenberg 
model, with a Curie like behavior at high temperatures 
and a broad maximum at a scale of the order of J, as 
the spins begin to form singlets. The main effect of dop- 
ing is to reduce the effective exchange. Experimentally 
the shift of the minima in the susceptibility is see n 131 i 132 , 
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FIG. 39: (Color online) The dynamical spin susceptibility at 
q = (tt, 7t) for few different doping levels and three differ- 
ent temperatures: superconducting state, normal state at the 
transition temperature and at room temperature. The pro- 
nounced peak is formed in SC state at 0.16t ~ 48meV and 
a broad peak in normal state is around 100 — 140 meV. Sus- 
ceptibility at normal temperature is much smaller and the 
peak moves to higher frequencies. The resonance is strongest 
at the optimally doped system. It disappears quickly in the 
overdoped side and somewhat more slowly in the underdoped 
side. Results are obtained with EDCA and NCA. 

but it occurs faster than in figure The effective spin- 
exchange will be reduced by the addition of a negative 
next-nearest neighbor hopping t' to the model. 

We now turn to the frequency dependence of the imag- 
inary part of the (jr, n) susceptibility probed in neutron 
scattering experiments. As shown in figure 1391 a pro- 
nounced peak in the (tt, tt) spin response at frequency 
0.16i in the optimally doped regime can be observed 
when entering the superconducting state. The position 
of the peak is temperature independent, but depends 
weakly on doping tracking the critical temperature. Our 
results are in qualitative agreement with experiment, for 
example the resonance energy scales with doping like 5T C 
and its position does not depend on temperature^. In 
addition we see a broader peak around 0.35 — 0.45£ ex- 
tending to very high frequencies of order of t ~ 300 meV 
which also gains some weight when entering the super- 
conducting state. 

Cluster methods coarse grain the momentum depen- 
dence. In the plaquette case, the coarse graining is done 
over 1/4 of the Brillouin zone centered at (tt, 7t) there- 
fore it is reasonable to compare our results with the q 
integrated susceptibility from Ref. 1 1331 where the two fea- 
tures, present in the mean field theory, 35 meV resonant 
peak as well as broader peak around 75meV extending 
up to 220 meV were observed. 

The exchange energy of the t-J model can be expressed 
as an integral of the spin susceptibility^ 

3 J f 

E xc = — d 2 qduj b(id)Im{x(q_, uj)](cos q x + cos q y ) 

(56) 
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FIG. 40: (Color online) The exchange energy versus temper- 
ature for few doping levels obtained by the NCA in EDCA. 
The lower panel is a blowup of the low temperature regime. 
The dotted lines correspond to metastable normal state be- 
low T c while the full lines continue into the superconducting 
state. 



Using this equation, one can elucidate the origin of the 
superconducting condensation energy, and the relative 
contribution of the different features of the spectral func- 
tion. 

Clearly, an important contribution to superconducting 
condensation energy arises from the incoherent features 
of the spin spectral function (around the frequency 0.4 — 
0.5i) rather than from the spin resonance. 

The exchange energy as a function of temperature is 
shown in Fig. 20] At high temperature, spins are disor- 
dered and the exchange energy is negligible. At temper- 
ature below J the singlets are formed and the exchange 
energy noticeably decreases, especially in the underdoped 
regime. At T c the exchange energy decreases further 
and gives far the largest contribution to the condensa- 
tion energy of the t-J model, as shown in Ref. H^. The 
exchange energy mechanism, observed in cluster DMFT 
study is thus in agreement with the strong coupling mag- 
netic mechanism for the superconductivity. 

The spin resonance ha s bee n viewed from two differ- 
ent perspectives (see Ref. 1 1291 and references therein) : i) 
starting from electronic quasiparticles and their residual 
interactions in a d-wave superconductor, residual inter- 
actions form a particle hole bound state with spin one, 
which is identified as the spin resonance, ii) Alterna- 
tively starting from a disordered quantum spin system, 
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FIG. 41: (Color online) left: Pseudoparticle spectral func- 
tions for the three most important pseudoparticles: ground 
states for N=4, N=3 and N=2 sectors. Right-top: Sketch of 
pseudoparticle threshold energies which can be interpreted as 
the effective many-body levels in normal and superconduct- 
ing state. Right-bottom: Pseudoparticle occupancies versus 
doping for most important pseudoparticles. The full lines cor- 
respond to the normal state while the dashed lines correspond 
to the superconducting state. 



one can identify the spin resonance as a massive spin one 
excitation, which becomes massless as one approaches the 
magnetically ordered phase. 

The cluster EDMFT equation [7] reconciles both points 
of view in a unified approach, since the equations for the 
spin susceptibility contain both the exchange interaction 
characteristic of the insulator J(q), as well as the quasi- 
particle contribution described by the spin cumulant M 
Eq.H 

The appearence of the spin resonance requires the dra- 
matic decrease of the anomalously large scattering rate in 
the normal state which is strongly reduced when the elec- 
trons condense to form d-wave pairs, avoiding criticality 
at low temperatures. The resonance, however, appears 
only in the superconducting state and is not present in 
the normal state. 



VII. PSEUDOPARTICLE INTERPRETATION, 
CONNECTION WITH OTHER MEAN FIELD 
THEORIES 

In this section we give interpretation of physical ob- 
scrvables in terms of pseudoparticle (cigenstates of the 
cluster) spectral functions. This is an alternative insight 
into a rich physics contained in the solution of cluster 
DMFT equations on a plaquette. 

Pseudoparticle creation and annihilation operators 
were introduced as mathematical entities representing 
the atomic eigenstates of the plaquette immersed in the 
Cluster DMFT medium. We have found that out of the 
large number (3 4 ) of pseudoparticles that we introduced, 
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very few of them are important for reproducing the low 
energy part of physical observables. For example, more 
than 95% of the one-particle spectral function at the 
Fermi level comes from a few convolutions of pseudopar- 
ticles in Eq. (|33|) . within NCA approach. This constraint 
is however not present for high energy part of the spec- 
tra such as Hubbard bands where contribution of most of 
pscudoparticles can be identified. The ground state and 
the low laying excitations are much more restricted and 
are a superposition of only a few atomic states. In the 
plaquette these important states are: 



r 4 = \n 


= 4,5 = 


0,K = 0) 


(57) 
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where N is the number of electrons in the cluster eigen- 
state, 5 and S z are the total spin and its z component of 
the cluster eigenstate, and K is momentum of the cluster 
eigenstate. 

Notice that although only few cluster eigenstates con- 
tribute to the ground state, the wave function is still 
highly nontrivial since it is a product state of an infi- 
nite number of states in the bath with the few atomic 
eigenstates of the impurity. This surprising result of re- 
striction to a few cluster eigenstates could be beneficial to 
devise useful approximations while extending C-DMFT 
to larger clusters in the future. In this paper we exploit 
this fact to give a simple interpretation of the different 
doping regimes of the t-J model. 

Fig. PHTleft shows the evolution of the three most im- 
portant pseudoparticle spectral functions from the un- 
derdoped to the overdoped regime. 

At small doping, the cluster is mostly occupied by the 
singlet state with one particle per site and zero momen- 
tum r 4 = \N = 4,5 = 0, K = 0), (half filled singlet). 
This pseudoparticle has the largest occupancy as shown 
in Fig. UTJ It describes a system locked in a short range 
singlet state as a consequence of the strong superex- 
change interaction. 

The cluster electron spectral function describes the 
process of addition and removal of an electron from the 
cluster at frequency to. Within NCA, it is constructed 
from the convolution of two pseudoparticles with differ- 
ent cluster occupation N and N + 1, or N — 1 , with 
the frequency restricted between zero and lo as described 
by Eq. (|3"3")) . The necessary condition for a peak of 
the one-particle spectral function at the Fermi level is 
that at least two pseudoparticle spectral functions share 
a common threshold and are strongly peaked at the same 
threshold. 

In the underdoped regime, the thresholds of all other 
pscudoparticles except r 4 are significantly shifted with 
reference to the half filled singlet, a pseudogap results in 
the one particle spectra in the underdoped regime. This 



gap in threshold energies severely limits the possible de- 
cay processes of the electron resulting in a low electronic 
scattering rate. This is the plaquette-impurity model of 
a few holes propagating in a sea of singlets. 

At large doping, i.e. in the overdoped regime, where the 
Kondo scale is dominant we obtain the standard DMFT 
description of a strongly correlated Fermi liquid. As is 
well known from the study of the Fermi liquid regime of 
the single impurity Anderson model, all pseudoparticles 
develop thresholds (X-ray singularity) at the same fre- 
quency which is related to the Kondo temperature of the 
problem. In our plaquette DMFT, all three important 
pscudoparticles (half-filled singlet, doublet with one hole 
per plaquette r3 CT K and singlet with two holes per plaque- 
tte T2) have a power law divergence at the same thresh- 
old frequency at zero temperature (Fig. I41[) which is a 
standard signature of the Kondo effect. Hence, the one 
particle spectral function begins to develop the Kondo- 
Suhl resonance at the Fermi level since the convolution 
between the doublet and half-filled singlet T 4 (or T2 
singlet) state is large at low frequency. The one-particle 
spectral function is peaked slightly above the Fermi level. 
Notice that while we cannot follow the formation of the 
Kondo resonance to very low temperatures due to the 
well known NCA pathologies, we can clearly see the on- 
set the Fermi liquid behavior in Fig. 1131 The overdoped 
regime characterized by the common threshold of pseu- 
doparticles is distinctively different from the underdoped 
regime, where only important state is the half-filled sin- 
glet r 4 and the r3 CT K doublet. The latter has a very little 
spectral weight in the region of the singlet peak. 

The Transition Region, Normal State: In the optimal 
doped regime, the Kondo scale and the superexchangc 
compete giving rise to a regime with very large scattering 
rate and consequently a small coherence scale. 

Surprisingly the evolution of the spectral function with 
doping is such that the optimally doped regime is approx- 
imately particle hole symmetric. As shown in Fig. 1411 
right-top the threshold of the N=2 cluster ground state 
and N=3 cluster ground state (doublet) merge first re- 
sulting in a Kondo-like contribution to the electron spec- 
tral function. This contribution is peaked above the 
Fermi level in a one band model below half-filling, in 
a Fermi liquid regime. The half-filled singlet however 
remains the lowest state in energy and still gives a sig- 
nificant contribution to the electron spectral function. 
The later contribution is peaked below the Fermi level 
and keeps a pseudogap-likc shape. Adding the two con- 
tributions to the electron spectral function restores the 
particle-hole symmetry in the density of states both in 
normal and superconducting state at optimal doping (see 
Fig. 113b and !25p . The approximate restoration is impor- 
tant, because it is known that cluster of impurities such 
as the two impurity Kondo modcl^, have a critical point 
only in the particle hole symmetric casein. 

Notice that the point of maximum scattering rate in 
Fig. [T7] coincides with the merging of the thresholds of 
the pseudoparticles (see Fig I4T1) . Around the same dop- 
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ing level, an approximate particle hole symmetry is re- 
stored in one-particle Green's function (see Fig. [5] and 
I13[) . Hence the term avoided cluster quantum multicrit- 
icality describes better the phenomena observed in this 
study since to reach the quantum critical point both the 
particle hole symmetry and the ratio of Kondo to RKKY 
coupling need to be varied. 

Transition into the superconducting state: The degen- 
eracy responsible for the strongly incoherent metal with 
large scattering rate at the Fermi level is lifted by the 
superconductivity avoiding the critical point. This dra- 
matic reduction of scattering rate in going from the nor- 
mal to the superconducting phase, depicted in figures fT7l 
and [18j highlights how anomalously incoherent the nor- 
mal state at optimal doping is, and how those anomalies 
are removed by superconductivity. 

This fact has also a natural interpretation in terms 
of pseudoparticles. Fig. HTHeft shows that both impor- 
tant singlet pseudoparticles (for T4 and T 2 ) develop a 
very sharp peak at the same threshold frequency and, 
at the same time, their occupancy increases (see Fig. WW 
bottom) upon condensation, indicating that electrons are 
locked into singlets with zero momentum. A gap opens 
between the singlets and doublets which gives the gap 
in the one-particle density of states. Because of this 
gap in the pseudoparticle thresholds, the large imaginary 
part of the electron self-energy does not persist in the 
superconducting state (see also Fig. [T7|) . Notice how- 
ever that in superconducting state the pseudoparticles 
are strongly mixed and the off-diagonal spectral func- 



tion ^4r 4 r 2 a l so develops a pole at the same threshold as 
Ar 4 and Ar 2 . The off-diagonal spectral function Ar 4 r 2 
describes the creation of a Cooper pair on the cluster 
Gr 4 r 2 = (0| a r 4 ( r ) a r 2 10) and therefore diverges at low 
temperature at the same threshold frequency. 

Since the density of states is composed of two almost 
equally important contributions, i.e., the convolution of 
the doublet with both singlets and T 2 ), the super- 
conducting gap is almost particle hole symmetric in the 
optimally doped regime with half-width of the order of 
O.lt. When the doping value is changed from its critical 
value, the asymmetry in the superconducting density of 
states appears. The magnitude of the asymmetry is the 
same as the asymmetry of the corresponding normal state 
spectra and comes from the fact that the occupancy and 
therefore importance of the T4 singlet exceeds the impor- 
tance of the r 2 singlet (see Fig. |41]-bottom) . 

Finally we comment on the role of the triplet pseu- 
doparticle. The spin susceptibility comes almost entirely 
from the convolution of the half-filled singlet with the 
half-filled triplet (T^ with ) . The later develops a peak 
at an energy 0.16i upon condensation which results in the 
resonance in the spin susceptibility shown in Fig. 1391 

It is interesting to derive the form of a low energy 
Hamiltonian involving the pseudoparticles in question. 
The conservation of charge, spin and cluster momentum 
considerably restricts the form of this Hamiltonian. If 
we assume it is of the Kondo form, it takes the following 
form 



-ff = ^ erapar + ^ e/ £ Qcj c q .CfcQ - + JiaJ, 4 ar 2 ^ e a!<T >CkQaCkQa' + h.c (62) 

r fcQa" kk'Qa' a 

+J2 ^2 "Lk^'k^^'k^'^Kt + h.c. + A)^(4flr - 1) 

kak'a'. K,K'e[(0,7r),(-7r,0)] T 

I 



where T runs over the relevant low energy pseudoparti- 
cles. e^cr' is an antisymmetric tensor and the Q runs 

over cluster momenta. Here c^q^ operators create an 
electrons in the bath with cluster momenta Q and spin 
a. The operators a r create a pseudoparticle on the clus- 
ter (see Eq. [T5]) . 

This Hamiltonian contains the competition of the par- 
ticle hole and particle particle channels for pairing with 
the baths of conduction electrons, and the approach to 
criticality is controlled by the variation of the on site 
energy er which should be identified with the pseudopar- 
ticle thresholds. It would be very intersting to inves- 
tigate this impu rity model with the tools and the per- 
spective of Ref. Il37l . It is clear that superconductiv- 
ity will add magnetic field like terms proportional to 
( a r 4 a r 2 ) Yshk'Qv'v W c feQ°-CfcQa-'- These terms should 



be strongly relevant, and moves the system away from 
criticality. 

Within CDMFT, the cluster of few sites (2 x 2 in our 
case) hybridizes with the Weiss field A, defined in Eq. ©. 
In single site DMFT, this effective medium drives the 
Mott transition. On the Bethe lattice within single site 
DMFT, it is proportional to the local Green's function 
A = t 2 G and therefore becomes gapped in the Mott insu- 
lating state while it remains gapless in the metallic phase. 
Hence due to the DMFT self-consistency condition this 
quantity shows a very strong doping dependence. 

Within cluster DMFT the effective medium is only 
weakly doping dependent and the evolution with doping 
is smooth (see Fig. |42|) in the doping range considered 
here. Moreover, this quantity shows very mild momen- 
tum dependence as opposed to strong momentum depen- 
dence of self-energy shown in Fig. [5) For examples, the 
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FIG. 42: (Color online) First panel: Imaginary parts of the 
cluster hybridization functions for various dopings in normal 
state at T = O.Oli using C-DMFT and CTQMC. Second 
panel: Real parts of the same hybridization functions in nor- 
mal state at T — O.Oli. Third panel: Imaginary parts of 
the same hybridization functions in superconducting state at 
T — 0.0054. Fourth panel: Real part of the same hybridiza- 
tion functions in superconducting state at T = 0.0054. The 
self-energies of the cluster show strong momentum depen- 
dence while hybridizations are only weakly momentum de- 
pendent. Furthermore, there seems to be no indication of any 
criticality in the hybridization functions such as the formation 
of a gap. 



(0, 0) and (%, w) Green's functions show almost no spec- 
tra at low frequency (are almost gapped) while the hy- 
bridization functions of these two orbitals are very similar 
to (0, 7r) hybridization function which contains most of 
the low frequency spectral weight. The mild and smooth 
doping dependence of hybridization functions leads us to 
believe that the proximity to quantum cluster criticality, 
which manifests itself in large scattering rate and van- 
ishing coherence scale, is driven by the impurity model 
itself rather than the self consistency condition. 

The picture here, is based entirely on a finite tempera- 
ture analysis, and is in the spirit of the DMFT approach, 
where we approach the strong correlation problem start- 
ing from high temperatures. 

It is important to continue the normal solutions of the 
plaquette DMFT equations to very low temperature to 
clarify the mathematical source of the criticality that we 
observe at higher temperatures. The critical point could 
occur exactly at T = 0, as proposed by Capone et. 
in the context of the two band Hubbard model with in- 
verted Hund rule exchange, and by us in ref^. The quan- 
tum critical point may exist in a impurity model with a 
fixed bath or might require the DMFT self consistency 
condition. Alternatively, there may be a finite second 
order endpoint of a first order line, as found in DMFT 
lattice models related to the two impurity models. No- 
tice also that power laws in an intermediate asymptotic 
regime, without an obvious underlying quantum critical 
impurity model have also been found in impurity model 
related to frustrated magnets^. 

Still, the precise nature of the low temperature normal 
state phase below T c is not essential for the validity of 
the CDMFT description. What matters is that at very 
high temperature T > J single site DMFT is a good de- 
scription of the system, and as we lower the temperature 
we find a broad region of temperatures where the pla- 
quette reference frame correctly captures the physics of 
the problem with its apparent criticality, even though at 
much lower temperatures a more non local description 
will be needed. It is even possible that the zero temper- 
ature solution of the DMFT equations do not exist, in 
which case, unlike the standard BCS theory where the 
superconductivity is viewed as an instability of a normal 
phase, we would have a superconducting state that exist 
without an underlying normal state. 

There is an important distinction between our views 
and those of local quantum critical scenario based on 
single site EDMFT scenario^. In the latter case, the 
locality of the quantum critical theory of the lattice is 
asserted to be reliable at T = 0, while the results of the 
EDMFT equations in two dimensions are known to be 
less reliable as temperature is raised^ . On the other 
hand, the results of plaquette DMFT are expected to 
become more accurate as the temperature is raised. 
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VIII. CONCLUSION AND DISCUSSION 

In this paper we developed and applied a plaquette Dy- 
namical Mean Field Theory to understand the nature of 
the superconductivity near the Mott transition. In rela- 
tion to earlier works, we focused of low but finite temper- 
atures to allow a comparison with the underlying normal 
phase. For this purpose advanced impurity solvers were 
brought to bear on the solution of the CDMFT equations. 

The idea of using a plaquette in a self consistent 
medium as a reference frame to reconstruct physical cor- 
relations functions on a lattice, while appealing, has sev- 
eral different implementations through different cluster 
schemes. Here, we stressed the numerous qualitative fea- 
tures which are common to all methods while pointing 
out the few significant quantitative discrepancies that we 
found among the different cluster methods in the course 
of our investigations. 

The low temperature landscape of strongly correlated 
electron systems can have many competing phases, for 
example conmensurate and incomensurate condensates 
of charge spin and current. A first step towards under- 
standing this landscape is to follow the evolution of well 
defined phases as a function of control parameters. In 
this paper we focused on the superconducting and nor- 
mal phase. Other phases and the competition with su- 
perconductivity can be studied with CDMFT techniques, 
as was done for example in Ref. [3?] for the conmensurate 
antifcrromagnetism. 

We find that the normal state in the mean field theory 
has two distinct regimes, which are naturally character- 
ized in terms of the regimes of the impurity model. At 
low doping, in the immediate proximity of the Mott in- 
sulating state, we have a realization of the RVB picture 
of holes propagating in a sea of spins with strong singlet 
correlations. In the impurity model language that corre- 
sponds to the RKKY phase of the two impurity model 
and its generalization to a plaquette. At high doping, 
we have a regime with well formed quasiparticlcs with a 
Fermi surface containing 1 — 6 electrons. In the impurity 
model language this corresponds to the Kondo regime of 
the one impurity model, and single site DMFT provides 
an adequate description of its properties. 

Plaquette CDMFT has three independent cluster self 
energies. For very large doping, only the local cluster 
self energy is non zero, indicating the validity of single 
site DMFT. As doping is reduced, S TW (iw) acquires a 
large real and imaginary part. This is controlled by the 
existence of a pole which approaches zero frequency at 
certain doping 81 (Sf ~ 0.1 in the t-J and 5f = in 
the Hubbard model). When combined with the cumu- 
lant periodization, this anomalous growth, gives rise to a 
topological transition associated to the formation of line 
of zeros in the Greens function (line of poles in the self 
energy) at zero temperature discussed in Ref. [43|. We 
call the doping at which the topological transition of the 
Fermi surface happens 8$. Notice however, that from a 
CDFMT perspective which focus on the finite tempera- 



ture description, an infinite self energy is not necessary, 
and all that is required to generate the pseudogap regime 
with its concomitant formation of Fermi arcs, is a self en- 
ergy which exceeds the bandwidth. 

We identified another critical doping, 6% associated 
with a maximum in the scattering rate of the third clus- 
ter degree of freedom Eott- This is an example of cluster 
quantum multicriticality. Namely, a mapping of a lattice 
model onto a quantum impurity model with a critical 
point. This critical point satisfies the following condi- 
tions: a) it requires a cluster of impurities for its exis- 
tence, hence it has no analogy in single site DMFT, and 
b) has at least two unstable (relevant) directions (for ex- 
ample the ratio Jkondo/ Jrkky and particle hole sym- 
metry breaking in the two impurity model Varma Jones 
critical point^). 

At a critical doping 8^ 1 there is an avoided critical 
point in the normal phase, which is near the doping level 
with the highest superconducting transition temperature. 
Since 8% > 5%, 5% may lie very close to 8$ if one adopts a 
periodization scheme along the lines of the cumulant pe- 
riodization, but this issue is left for future studies since 
it strongly depends on the periodization scheme. Look- 
ing at the scattering rate and coherence temperature in 
the normal state solution of CDMFT equations, we thus 
identified a critical doping 8%, which could be related to 
the hidden quantum critical point which has been hy- 
pothesized by many authors based on a large body of 
experimental dat a 135 i 142 . 

We have not analyzed the properties of the CDMFT 
quantum impurity model describing the normal state at 
zero temperature. It is even possible that the normal 
state solution simply does not exist at T = 0. These 
points are largely academic from the point of view of 
the finite temperature physics which we want to describe 
with CDMFT. The manifestations of the possible quan- 
tum criticality are rapidly removed by the onset of super- 
conductivity. The electronic lifetime on the Fermi surface 
controlled by im£(0, it) is dramatically reduced as the 
system goes superconducting. 

One then arrives at a superconducting state, which 
inherits the normal state gap, largely caused by S„, 
but with coherent nodal quasiparticles characterized by 
a weakly doping dependent velocity perpendicular to the 
Fermi surface. The velocity along the Fermi surface 
va has a dome like shape and decreases in the under- 
doped regime providing further support to the two gap 

picture of the superconducting state of the underdoped 
cuprates^^^^. 

The superconducting phase is stabilized by the gain 
of superexchange energy, namely improved spin-spin sin- 
glet correlations. We resolved the frequency dependence 
of the anomalous Greens function and the anomalous self 
energy, and found them to have a structure very different 
from conventional phonon mediated superconductivity in 
the Migdal Eliashbcrg theory Since the superconducting 
state restores coherence, long lived sharp excitations, Bo- 
golubov quasiparticles and a sharp spin mode which re- 
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scmbles the neutron " 40 meV resonance" , emerge below 
T c . 

We extracted different observables such as tunneling 
density of states, optical conductivity, optical mass and 
plasma frequency, integrated optical spectral weight su- 
perfluid stiffness and spin susceptibility which compare 
well at a qualitative level with experimental data on cop- 
per oxide materials. 

We have shown that at 5^ (which occurs very near the 
maximum in T c ) the coherence energy vanishes and the 
scattering rate is maximal. At this doping, an approx- 
imate particlc-holc symmetry in the one electron spec- 
tra is recovered, and approximate power laws in physical 
quantities (a oc cj~ 2 / 3 ) emerge in an intermediate fre- 
quency range. 

Upon periodization, the large value of the non local 
self energies turn Fermi surface into Fermi arcs^ 3 -, and we 
studied the evolution of the fermi arcs with temperature. 
We showed that within C-DMFT fermi arcs shrink with 
decreasing temperature. 

Our solution of the CDMFT on a plaquette, has many 
similarities with the earlier studies of Anderson's Res- 
onating Valence Bond theory of high temperature super- 
conductivity in the slave boson mean field theory for- 
mulation. This approach, correctly predicted the d-wave 
symmetry of the superconducting order parameter and 
the presence of a pseudogap with the same symmetry 
well above Tj&. 

The similarity between the results of the CDMFT and 
slave boson approaches is not accidental. Both methods 
are mean field techniques based on order parameters that 
can be defined within a plaquette, and capture the effects 
of the proximity to a Mott insulating state using a small 
set of short range degrees of freedom. 

Compared with slave boson mean field theory, CDMFT 
has additional flexibility both in the frequency of the 
one electron spectral function as well as in its momen- 
tum dependence. One crucial difference, is a much more 
pronounced momentum space differentiation with very 
different electronic properties at the nodes and at the 
antinodes. This anisotropy, with the concomitant ex- 
istence of two energy scales in the superconductor, re- 
solves the earlier problems of the RVB theory related to 
the doping dependence of the linear term of the pene- 
tration depths. The need for the introduction of more 
anisostropy in the microscopic theory had been antici- 
pated by experiments, and by the phenomenological anal- 
ysis of Ioffe and Milli a 144 ' 145 ' 146 . Recent phenomenolog- 
ical model o 82 i 83 have also generated a more pronounced 
momentum dependence of the one particle spectra, and 
incorporated in their approach a va that decreases with 
decreasing doping. The main differences between these 
phenomenological approaches and the more microscopic 
C-DMFT reside in the location of the lines of zeros of 
the Green function. While in Refs. HHH! the lines of zero 
lie on the the umklapp surface, in C-DMFT the lines of 
zeros are dynamical entities with a location that evolve 
with doping. 



CDMFT is an extension of single site DMFT an ap- 
proach that has been very successful in describing many 
aspects of the finite temperature Mott transition. Using 
a single impurity in a medium, this method, has been 
able to describe several regimes near a Mott transition. 
A Fermi liquid regime, at small U and temperature, a 
bad metal at temperature larger than a characteristic 
temperature T co h(U), a Mott insulator at large U and 
low temperature, and a bad insulator when the temper- 
ature is high enough that the Hubbard bands begin to 
merge. 

By construction, single site DMFT assigns the same 
self energy to the electronic states on the whole fermi 
surface. Hence at a given energy and temperature, cither 
all the states at all k points are coherent, or they are all 
incoherent. This is not a good description of the high 
temperature superconductors, which therefore cannot be 
described with single site DMFT. 

On the other hand, CDMFT allows the states in the 
nodal region to be coherent quasiparticlcs while at the 
same time the states in the antinodal region are highly 
incoherent and have a pseudogap, i.e., T co h(nodal) <C 
T co h(antinodal) . The self energy in the nodal region 
could be compared to a single site DMFT in the Fermi 
liquid regime with U < U C 2 and T c < T co h{U) while the 
antinodal self energy is more of a single site DMFT in the 
bad insulator regime U > U C 2 and T c > T co h{U). Pla- 
quette DMFT, offers a mean field picture of the lattice 
problem, whereby the different cluster self energies and 
cumulants describe different regions of momentum space 
with distinct physical properties. A nodal region which 
is closest to a Fermi liquid, an antinodal region which ex- 
hibits a pseudogap, and an intermediate region between 
the two, described by the (0, 7r) self energy, which exhibits 
the maximum scattering rate at criticality. 

This qualitative picture is only a crude caricature of 
the full CDMFT solution, but it is a useful qualitative 
guide to understand how the Fermi arcs originate from 
the proximity to the Mott insulator, and above all, high- 
lights why single site DMFT is inadequate in this situa- 
tion. 

The objective of this work was to advance our under- 
standing of the t- J and Hubbard model as a " bare bones" 
model of the density driven Mott transition. Important 
open problem is to incorporate and understand how other 
effects, such as the effects of more realistic band struc- 
ture in the multi band model, the disorder and the elec- 
tron phonon interaction, which play an important role in 
cuprates, can effect the solution of the model, 

We presented a qualitative comparison with several ex- 
periments in the copper oxide based materials, and given 
the limitations of the model and of the methodology used, 
this comparison is very encouraging and warrant future 
studies including more accurate modeling and further 
methodological improvements. 

Future studies should include a realistic band struc- 
ture of the copper-oxygen planes and additional Coulomb 
terms beyond the local Hubbard U, that can be accomo- 
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dated on the plaquette. In addition to d x 2_ y 2 copper 
band, it would be desirable to include another copper 
band, namely, d z 2 band which is coupled to apical oxy- 
gen. Although the latter band is filled in the band struc- 
ture calculation, it comes close to the Fermi level. 

Another important direction is to better momentum 
resolve one particle and two particle quantities. The 
latter will require advances in the analytic continuation 
techniques of QMC data, as well as a better understand- 
ing of how to convert cluster quantities into lattice ob- 
scrvables in C-DMFT. Furthermore within a cluster size, 
it is important to implement an optimal choice of orbitals 
in CDMFT describing different momentum patches in the 
Brillouin zone. Functional approache o 23 i 147 as well as 
CDMFT inspire d, m odeling of experimental data along 
the lines of Ref. Il48l can provide useful guidance in this 
direction. 

Mean field approaches clearly separates the short dis- 
tance effects contained in the theory, from long distance 
effects which will require the introduction of fluctuations 
due to vortices and pair fluctuations. The T c vs 8 line 
in CDMFT should be interpreted as being close to the 
Nerst line in the cuprate phase diagram^. On the other 
hand, the true superconducting critical temperature line, 



is strongly reduced relative to the CDMFT on the un- 
dcrdoped side of the phase diagram to the effects of long 
wavelength fluctuations of the order parameter, which re- 
quire long wave lengt h field theoretical techniques along 
the lines of Ref. Il50l 

Finally other inhomogeneous phases, such as stripes, 
charge, bond, pair density waves, and other broken sym- 
metries can appear as secondary instabilities, and can be 
studied with our methods by inserting a relatively local 
(restricted to a plaquette), but site dependent self ener- 
gies into the CDMFT functional. 
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